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We are stardust 
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Abstract 

This paper reviews the formation of carbon dust in the outflow from stars and 
the subsequent evolution of this socalled stardust in the interstellar medium. The 
chemical and physical processes that play a role in carbon stardust formation are 
very similar to those occuring in sooting flames. Based upon extensive laboratory 
studies of the latter, the structure and physical and chemical properties of carbon 
soot are reviewed and possible chemical pathways towards carbon stardust are 
discussed. Grain-grain collisions behind strong interstellar shocks provide the .high 
Dressures required to transform graphite and amorphous carbpn grains m *c diamond. 

^fprocessls examined and the properties of shock-synthesized diamonds are 

reviewed. Finally, the interrelationship between carbon stardust and carbonaceous 
meteorites is briefly discussed. 


I Introduction 


Interstellar dust is an important constituent of the interstellar medium. It 
absorbs efficiently visible and ultraviolet photons, which it reemits at R 
wavelenaths and thus it regulates the spectral appearance of many dust-shrouded 
objects. Tt regulates the energy balance of the gas through the ejection of energetic 
electrons and gas-grain collisions. By attenuating the ambient interstellar UV 
radiation field, it allows (indirectly) the formation of quite complex molecules. 
Undoubtedly, the most abundant molecule, Ho, is formed directly on grain surfaces. 

The abundance of other molecules may also be influenced through gram-catalyzed 
reactions. Grains may also have played an important role in the cosmic history of the 
biogenic elements (Wood and Chang 1985). In fact, carbonaceous grains are thought to 
form an important component of interstellar dust (1). 


In recent years, our understanding of interstellar carbonaceous dust has 
improved dramatically. Partly, this is a result of the opening up of the infrared 
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window Infrared spectroscopy is a major tool for remote probing of the composition 
of interstellar dust (2). Indeed, most dust components thought to be present in the 
interstellar dust have been identified through infrared spectroscopy. The development 
of new techniques that allowed detailed analysis of microscopic amounts of stardust 
preserved in carbonaceous meteorites has also been instrumental in increasing our 
knowledge of interstellar dust (3). y 


The picture that arises from this new wealth of data is one where interstellar 
dust consists of many diverse components formed under a variety of conditions in 
many classes of objects. This varied background of interstellar dust makes the 
quotation from Woodstock particularly appropriate. A multitude of complex processes 
have been suggested to play a role in circumstellar dust formation and its 
modification in the interstellar medium. In view of the rapid developments, it has not 
vet been possible to weed out the unimportant processes and to concentrate on the 
tew, more important ones. Lacking general consensus in the field, this review is an 
attempt to summarize my own views on these problems. Some of these seem to eniov 
considerable following, others are more controversial. However, given the current 
large flux of ideas, any review has to be somewhat personal and dated anyway. 


This review is organized as follows: 

I Introduction 

II Composition and Evolution of Interstellar Dust 

III The Structure of Carbonaceous Materials 

IV Carbon Injection into the Interstellar Medium 

A) Sources 

B) Carbon budget of the galaxy 

C) The 13 C budget 

V The Dust Condensation Zone 

A) The photosphere 
Bj Structure ot the outflow 

Cj Physical conditions in the dust condensation zone 

D) Chemical composition 

VI Chemical Pathways to Carbon Stardust 

A) Thermodynamics and carbon condensation 
Bj Soot formation 

CJ PAHs and the formation of carbon stardust 

D) Non-LTE and the formation of carbon stardust 

|) Carbon stardust formation around WC 8-1 0 and R Cr B stars 

F) Carbon stardust formation in supernova ejecta 

G) Summary 

VII Grains and Shocks in the Interstellar Medium 

A) Grain-grain collisions 
Bj Diamond metamorphism 
C) Diamonds in the sky 

VIII Interstellar Carbon Dust and Carbonaceous Meteorites 

IX Summary 


Composition and Evolution of Interstellar Dust 
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Interstellar dust consists of many different components, including amorphous 
silicates, graphite, Polycyclic Aromatic Hydrocarbons (hereafter PAHs), amorphous 
carbon, icy grain mantles, and organic refractory grain mantles. Stoncarbide and 
maqnesium sulfide, which have been identified in the outflow from carbon- rich giants 
and planetary nebulae might also be present in the interstellar medium. Our present 
astronomical knowledge on the composition of interstellar dust is summarized in 
table 1 For each component, this table gives the birth site(s), the faction of the 
elemental abundance that it locks up, and the dust volume relative to that of the 
silicates. The spectral signatures that have been used to identify each dust 
component are also indicated. For a discussion of these the reader is referred to a 
recent review by us (1) as well as to earlier reviews (4,5,6). 

One of the most interesting developments in recent years within the field of 
interstellar dust has been the realization that large polycyclic aromatic molecules 
form an ubiquitous component of the interstellar dust. A variety of objects show 
prominent emission features at 3.3, 6.2, 7.7, and 11 .3 pm (5,7,8), which are 
characteristic of PAHs either as free-flying molecules or in the form of soot 
Darticles (9 10,1 1). Figure 1 compares the IR spectrum observed towards the HI I 
region in Orion with the Raman spectrum of auto soot (11). For a highly disordered 
material such as soot the infrared active and Raman active vibrational modes are 
similar in number and frequency. It is expected that a collection of free ’ .^y* n 9 ^ 
molecules will have a similar IR spectrum. Although it is disti^^^ 

spectrally between a collection of large PAH molecules and solid soot particles, an 
analysis of the excitation mechanism shows that the carrier of these emission 
features contains only about 20 C atoms (10,1 II ,12). Depending on the FUV ^bsorp !on 
properties of the particular PAHs responsible for the IR 

total available elemental carbon is locked up in this component (cf., table 1)- Since in 
many objects the IR emission originates in material recently ejected fro ™ ® ®tar, 

PAHs probably are the extension of the grain size distribution of circumstellar carbon 
gratos into the molecular domain and as such may shed light on the deta. ed formation 
mechanism of carbon stardust. This will be discussed in some detail in § VI. 


About half of the volume of the interstellar dust is made up of silicate grains 
(c f table 1 ). The other half has to be in the form of a largely carbon-bearing dust 
component/either graphite or amorphous carbon stardust, or organic refractory grain 
mantles accreted on stardust in the interstellar medium. The latter component 
consists of a mixed organic polymer containing CH 3 , CHo, C=0, and possibly OH groups, 
which is formed by UV phtolysis and (transient) heating of an icy (eg., H 2 0, NH 3 , 

CHoOH and CO) mixture. Current observations indicate that amorphous carbon 
particles seem to make up only a relatively minor fraction of the intersteHar dust 
Fc f table 1), although this value is somewhat uncertain. From table 1 , one might get 
the impression that the carbon dust component is dominated by organic refractory 
arain mantles. This is somewhat misleading. First, large graphite particles (>1000A ) 
do not have readily discernable absorption features and they are, thus, not accounted 
for in table 1 . Second, the relative volume of the organic refractory gram mantles 
which has presently only been detected towards the galactic center, is uncertain by a 
factor of about two due to the uncertainty in the precise identification of the 
absorbing molecules. Moreover, the relative volume of this dust component in the 
local interstellar medium may be much less. In fact, observations towards n ©3rby 
stars yield an upper limit on the relative volume of organic gram mantles which is 
about half that quoted in table 1 . In view of these uncertainties, the distribution of 
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FLUX (1<T 17 WcnT 2 pm) 


MICRONS, pm 

10 8.3 7.14 6.25 5.55 



Figure 1 • Comparison of the 5 to 8 pm Raman spectrum of auto soot with the emission from the Orion Hll 
region. 



Table 1 : Interstellar Dust Components 


COMPONENT 

BIRTH SITE 

ELEMENTAL 

ABUNDANCE 3 

REL VOL b 

SPECTRAL SIGNATURE 

Silicates 

O-rich giants 
& Novae 

100% Si 
20% 0 
(Mg & Fe ?) 

1 

10, 20pm features 

Graphite 

C-rich giants ? 

>25% C 

>0.25 

2200A bump* 

PAHs 

C-rich PN 

1%C 

0.01 

3.3, 6.2, 7.7, 11.3pm 
emission features 

Amorphous 

Carbon 

C-rich giants, 
Novae, R Cr B 
& WC 8-10 

5-1 0%C 

=0.1 

7.6pm absorption feature ? 

Icy grairi 
mantles* 3 

molecular 

clouds 

up to 40 % 
CandO 

up to 3 

3.1, 4.6, 6.0 and 6.85pm 
absorption features 

Organic 
refractory 
grain mantles e 

Interstellar 

medium 

24% C 
6% O 

0.8 

3.4, 6.0 pm absorption 
features 

SiC 

C-rich giants 
& Novae 

_c 

_c 

1 1 .4pm emission feature 

MgS 

C-rich Giants 
& PN 

_c 

_c 

30 pm emission feature 


a) Percentage of cosmic abundance of element locked up in interstellar dust component. 

b) Volume of dust component relative to that of silicates. 

c) This dust component has not (yet) been detected in the interstellar medium. 

d) Only observed in dense molecular clouds. 

e) Only observed towards the galactic center. 

f) Only graphite grains of 200A contribute to this bump. 


the carbon in the various carbon-bearing dust components in the interstellar 
medium should be considered approximate at this point. 


Detailed models for the dust in the interstellar medium based on both pre mises 
- graphite and organic refractory grain mantles - have been developed j. 1 f). ^ smg 

some combination of grain components, and a grain size (and shape) distribution, 
these models fit the visible and UV interstellar extinction curve the linear and 
circular polarization, the visual extinction per hydrogen atom and, to some extent, the 
observed elemental depletions. Such an exercise is useful for the derivation of 
important quantities such as the grain size distribution within the framework of one 
particular model. However, these fits cannot be used to discriminate between 
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EVOLUTION OF INTERSTELLAR DUST 


STARDUST : SILICATES, GRAPHITE, AMORPHOUS CARBON, PAHs, SiC, MgS 

FORMATION PROCESSES IN CIRCUMSTELLAR SHELLS: NUCLEATION, CONDENSATION, COAGULATION 



DUST IN THE INTERSTELLAR MEDIUM: STARDUST, ICY GRAIN MANTLES ORGANIC 
REFRACTORY MANTLES, DIAMONDS 

FORMATI ON PROCESSES IN MOLECULAR CLOUDS : ACCRETION, REACTION UV PHOTOLYSIS 
TRANSIENT HEATING 

DUST MODIFICATI ON PROCESSES : (DIAMOND) METAMORPHISM AND SHATTERING IN SHOCKS 


Figure 2 - A schematic representation of the evolution of interstellar dust, summarizing the processes that 
play a role in the formation, destruction and modification of interstellar dust. Stardust is formed in the 
outflow from stars in the late stages of their evolution and ejected into the interstellar medium. 

Interstellar medium dust is formed inside dense clouds. Dust is destroyed by strong shock waves and by 
incorporation into newly formed stars. Grain-grain collisions in shock waves can also transform graphite 
(or amorphous carbon) into diamonds. See text for details. 
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different qrain models. The visible and UV interstellar extinction curves are notably 
insensitive to the exact grain composition because the amorphous structure of 
interstellar grains as well as the presence of a size distribution tend to smear out 
any identifiable structure in the ejectronic absorption properties ; of the jaten^s (1). 

It is because of this that widely disparate models for the interstellar dust all yield 
reasonably good fits to the observed interstellar extinction curve. The presence of 
free parameters in the models such as the shape pr surface structure of the grains, 
whicn have not yet been considered in such exercises, reinforces this point. 

The evolution of interstellar dust is determined by many complex, interplaying 
Drocesses whose details are not understood very well. Figure 2 gives a schematic 
overview of the evolution of interstellar dust. Spme dust components -silicates, 
qraphite, and amorphous carbon - are formed in the outflow from stars (i.e., stardust) 
bv nucleation, condensation and coagulation processes and their composition and 
sizes reflect the chemical composition and physical conditions in the condensation 
zones In contrast, other dust components -icy and organic refractory grain mantles - 
are formed by accretion, reaction, UV photolysis and transient heating processes in 
the interstellar medium (i.e., ISM dust). Prolonged processing of such grain mantles 
may eventually lead to the formation of a highly disordered compound consisting 
mainly of carbon arranged in an aromatic network. This is akin to pyrolysis of 
hydrocarbons and will ultimately lead to an amorphous carbon dust component. Dust is 
also processed in the interstellar medium by strong supernova shpck waves. High 
velocity ion-grain collisions behind a shock will lead to sputtering, 'OJ-'ir p ’^ t f t ff’ rc 
and radiation damage. These will also lead to an amorpnization of the surface layers. 
Hiqh velocity grain-grain collisions send strong pressure waves into each grain which 
can lead to vaporization, shock metamorphism (eg., graphite -> diamond), and 
shattering. Since each of these processes will lead to a different grain material, the 
composition of the interstellar dust in a particular region will reflect the relative 
importance of these processes there. The composition of interstellar dust is 
therefore in some sense a misnomer since the dust composition is expected to vary 
with phase (eg., diffuse ISM versus dense cloud) as well as galactic radius (eg., local 
versus galactic center). A global outline of the physical processes playing a role in 
the formation and modification of interstellar dust has been given elsewhere (2) Here 
we will concentrate on the formation of carbon stardust and its subsequent evolution 
in the interstellar medium. 


Ill The Structure of Carbonaceous Materials 


Carbon-rich materials show a wide range of behavior due to the diversity in 
chemical valence structure and bonding possible with carbon. These are re ^® r '‘®d to as 
sn S p2 and sp^ hybridization, where the notation comes from the original carbon 
atomic orbitals involved in forming the molecular orbitals. In sp hybridization, the 
bond geometry is linear and two unequivalent bonds are formed, one 
single. An example of a simple organic compound with this structure is acetylene, 
H_C;=Q-H If only carbon is involved, an acetylene polymer (eg., polyynes) is formed, 
with a linear structure (cf., fig. 3). Note that polyacetylene is a polymer , denved from 
acetylene by opening up of the triple bond and has an alternating double bond 
structure (eg., sp^ bonding). In sp“ hybridization, the bond geometry is trigonal and all 
the bonds lie in a plane. In general only two of the three bonds are equiva ent and two 
single bonds and one double bond characterize this configuration. A simple organic 
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THE VARIETY OF CARBON SOLIDS 



sp HYBRID ACETYLENE POLYMER 




GRAPHITE 




sp 3 HYBRID 


DIAMOND 


E i gure 3 - A schematic representation of the chemical valence structures and bonding arrangements 
possible with carbon. Left: the orientation of the electron clouds in the different hybridizations. Right: the 
corresponding carbon solids. Solid lines indicate actual bonds, while dashed lines are intended to guide the 
eyes. Note that sp, sp^, and sp d hybridization lead to a linear, planar and three dimensional crystalline 
structure, respectively (see text for details). 
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compound with this structure is ethene, H 2 C=CH 2 . However, if only C ' s J nv0 !^®^ and 
there are enough atoms to form six-membered rinas, all the bonds can I 
pauivalent within the carbon framework. This results in an aromatic network as in 
araohite (cf fig 3). Hydrogen can bond to the external single bonds of the sd carbon 
atoms at the peripnerv of tne aromatic structure. PAHs are molecular examples of 
such structures. In spr* hybridization the bond geometry is tetrahedral and ^f a y 
carbon is involved a diamond structure is produced (cf., fig. 3). Examples of simple 
molecules with this structure are methane (CH 4 ) and ethane (H 3 0-OH 3 ). 

Any realistic material will have some impurities, defects and dislocations and 
pure qraphite, diamond or acetylene polymer materials do not exist >n nature Jn 
nartimlar a laroe variety of caroonaceous materials contain a mixture of C atoms 
S sp^and s£ bonding Thlse materials are all loosely called "amorphous - carbon" ,n 
the astroohvsical literature, although their physical properties can vary considerably. 
Hydrogenated amorphous carbon films are one example of such a |j° nding 

structure and soot is another. Both materials contain arge ara ^ at . 1 ^ dam . a,ns 
intersDersed with tetrahedral C atoms which cross link aromatic planes. In 
hydrogenated amorphous carbon these tetrahedral C atoms arerando ™' y S' bu n f d 
between the aromatic C atoms. In soot they are more concentrated at the edges of 
aromatic planes (cf. § VIB). In both cases a large concentration of hydrogen atoms 
(H/C=0 1) saturates’the tetrahedrally carbon and is bonded to the of the 
aromatic domains. This difference in amorphous character is ess p at ' a i ly 
between a "truly" (random network) amorphous material and a polycrystallme 

material (cf., 1). 

The physical and chemical properties of (hydrocarbons are 'argeiy determined 
bv these bondinq arrangements. For example, due to the conjugation of n electrons 
(ea sharinq of rc electrons between different C atoms), graphite shows a metalic 
behavior within the basal (aromatic) plane (eg., graphite is a semi-metal). Its optical 
DroDerties from the FUV throughout the IR are therefore dominated by interband and 
mtraband transitions of these "free" electrons (1 8). When the aromatic structure is 
broken up by the presence of hydrogen atoms or tetrahedrally bonded C atoms, ^ ’ 1 
amorphous carbon, then the electrons will be more localized. As a result, amorphous 
carbon is a semi-conductor with a band gap up to 5 eV depending on the size of the 
laroest aromatic domains. While its UV properties are still dominated by interband 
electronic transitions, its IR properties reflect now vibrations of the carbon (and 
hvdroaen) skeleton (1). It should be emphasized that many of the properties of solid 
carbon, in particular the electronic structure, are dependent on the disorder and not 
the disordering agent present. 


IV Carbon Injection into the Interstellar Medium 


A Sources 

There are many stellar sources that contribute to the cabon in the interstellar 
medium including caroon-rich red giants, carbon-rich Wolf Rayet stars, and 
suoernovae. The relative importance of these different types of objects is not wel * 
known Even estimates of the relative importance of hign mass stars versus low mass 
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Table 2: Stellar Mass-Loss Rates [M @ yr’ 1 kpc* 2 ] 


Object 

M H a 

M c b 

M d c 

Observational evidence for dust 

C-rich 

giants 

5(-4) 

5(-6) 

2(-6) 

1 ) featureless continuum: amorphous carbon 

2) 1 1 ,3pm feature: SiC 

3) 3.3,6.2,7.7,11.3pm features (in PN 
phase): PAHs 

O-rich 

giants 

5(-4) 

2(-6) 

- 

1)10 and 20pm features: silicates 

Novae 

6(-6) 

4(-8) 

1(-8) 

1) 10pm feature: silicates 

2) black body (marbles ?) 

3) black body + 1 1pm feature (SiC ?) 

OB stars 

1-3(-4) 

1(-6) 

- 

1 ) no dust 

SN: Type 1 
Type II 

2(-4) 

3(-7) 

2(-6) 

1(-6) 

1) no conclusive evidence for dust 

WR stars 0 * 

- 

4(-6) 

- 

1 ) no dust 

WC 8-10 


2(-6) 

5(-8) 

1) featureless continuum + 7.7pm 
feature: amorphous carbon (+PAHs ?) 

TOTAL 

1.5(-3) 

1.6(-5) 

3(-6) 



notes: 

a) Hydrogen injection rate. 

b) Carbon injection rate. 

c) Carbon stardust injection rate. 

d) Excluding WC 8-10 stars. 


stars vary continously due to new developments in, for example, the determination of 
mass-loss rates (eg., IR studies) or nucleosynthetic reaction rates (eg., 12 C (cfvW 
). Thus, while 10 years ago high mass stars (M>10 M J where thought to be mainlv 
responsible for the carbon in tne interstellar medium (?9), in more recent studies ^ 

th a t tho non ^ ian * s saem to dominate (20,21). However, there are some indications 
^i he , pen ?^ um IS f wm 9 |n 9 back into the direction of the massive stars In the 

3 nwSST ° thlS . sectian v astimate ? of the carbon mass-loss rates of different types 
of objects are reviewed. These estimates are summarized in table 2. P 

Bed Giants : The mass injection rate of gas in the interstellar medium is dominated hv 


* u birth ratesAss^ dwarf has^ 

lost about 1 M® during the red giant phase, we find from the birth rate of such obi 


objects 
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(=1 V r 1 - 22) a total mass injection rate of about 1 yr 1 , corresponding to about 
lO' 2 ^ kpc* 2 yr* 4 (23). This is supported by IR and CO studies of the dust ^nd ga^ 
mass-loss rate from these objects, which yield values between 5 and 9x10 M 0 yr 
k DC -2 (21 24-27). However, estimates of the mass injection rates due to Cwicn ana 
(j-ricn qiants separately range from equal amounts by the two types of objects 
(21 26) to O-rich giants dominating by a factor ten (24). The giants with a high 
mass-loss rate (>10* 8 M 0 yr 1 ) dominate the total gas return to the galaxy and the 
above mentioned discrepancy centers on differences in the total nu, ^ e . r h ° . 8 V, C ^ 

C-rich objects estimated to be present in the galpy. It should be noted that an 
studies require large correction factors for sample incompleteness and the results 
are thus inherently uncertain. We will adopt a return rate by C-riph [Qiants of 5x10 
yr 1 kpc* 2 , which is about the maximum estimated (25 27). Estimate^ on the 
refurn rate by O-rich giants agree on approximately 5x1 O' 4 M @ yr kpc (24 27)^ 
Using these gas return rates we have estimated the C return rates (cf^, table 2) by 
assuming an elemental C abundance of 5x1 O' 3 and 10. by mass for .^' n .9j 1 k" r c 

qiants respectively. The former corresponds to solar abundances, while the latter 
reflects the P C enhancement due to dredge up of freshly synthesized material during 
the AGB phase (cf., § VA). 

OB stars: Estimates of the total mass-loss rate of OB stars range from 1 1 to 3*10 4 
M vP kpc^ (28 29). With a solar elemental C abundance this translates into about 

1(FM @ yr 1 kpc* 2 . 

Novae : Nova events are generally interpreted as a thermonuclear runaway in ' an 
nwaTafe of and an^|cfo^ ^ 

'31-34). Assuming a C abundance of 7x10^ by mass leads to 4x10 8 M @ yr 1 kpc 
! or the C ejection rate (cf., table 2). 

Bunernovae: There is little direct observational information on the carbonjnjection ^ 
rateduetosupernovae (35) and their contribution is largeKr guessed at. There ar^ two 
types of SN, whose only connection is the large amount of energy released (~10 erg) 
and the explosive nucleosynthesis of heavy elements. Type l[ is general y believed to 
result from the gravitational collapse in massive stars (M > 8 M e ). Type I s are 
thought to be due to a thermonuclear instability of an accreting white dwarf by the 
accumulation of a critical mass (36). The total supernova ra t e Ti a s bee ne s t i m a ted to 
be about 2x1 0‘ 8 kpc*^ yr* 1 with about equal rates for type I and II (37). The typical 
type II supernova has a main sequence mass of 25 M e of which about 2 M 0 remains 
behind in a neutron star (36), this yields a total H mass injection rate of about 
2x1 0* 4 M vr* 1 kpc* 2 . Such a supernova ejects about 4 M @ in the form of freshly 
synthesized heavy elements, of which 0.25 M 0 is in the form qj T h |. C T , u i l | ctl0n 

rate from type II supernovae is then about 2x10 M 0 yr kpc (cf., table 2). This 
number is very uncertain since it depends critically on the somewhat uncertain C 
oc reaction rate and its effect on the presupernova convective shell structure 

36,38). 


t 


Type I supernovae will contribute much less to the galactic enrichment of 
intermediate-mass elements such as C. Calculations typically yield a C/Fe mass 
fraction ratio of about 3x1 0* 2 of the solar ratio (39). Since about 0.5 M @ of Fe is 
required to explain the lightcurves of type I supernovae (36), this implies the 
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production of about 2.5x1 O’ 2 M 0 of C per type I SN explosion. The resulting C ejection 
rate for type I supernovae is then about 3x1 0 7 M @ yr f kpc' 2 . 

Stars : Popu latio n I Wolf Rayet stars are the descendants of massive 
(M >25-40 M 0 ) OB stars where substantial mass-loss during the H burning phase as 
well as mixing processes have "uncovered" enriched stellar core material on the 
surfaces of these objects (40). The WC stars contain He, C and O and no evidence for H 
or hJ, while the WN have greatly enhanced He apd N and no H (41). The total mass-loss 
rate of WC stars is estimated to be about 3x1 O' 5 M & y r 1 kpc’ 2 where about 25% is 
due to WC 8-10, which are associated with dust formation (42). Observations of the 
C/He abundance ratio in these stars lie in the range 0.1 to 2 by mass (41 43) while 
theoretical calculations yield about 0.5 (44). Assuming^ a conservative value (0 2) for 
this ratio yields a total C injection rate of 6x1 O’ 6 M. vr* knr-2 f or a || wq s ^ ars an ^j 


B) Carbon budget of the galaxy 

From this discussion it is obvious that, although the mass return rate is dominated by 
red giants, many different object contribute substantially to the carbon budget of the 

' . uncertainty in these estimates should, however, once more be emphasized 
While previously massive stars where thought to contribute mainly during the 
supernova phase (19), the upward revision in the 12 C(a,y) 16 0 reation rate makes this 
a less important source of C. Nevertheless, recent observations show that massive 
stars still contribute substantially albeit during the Wolf-Rayet phase. Essentially 
during this phase they mix their helium-burning products to the surface and eject 
them in a strong wind before further nucleosynthetic burning can take place. 


Elemental abundance studies of external galaxies may provide additional constraints 
on the type of objects that produce the C enrichment of the galaxy (45). For example 
the low C/O ratios observed in the Large and Small Magelanic Clouds show that since 
O is mainly synthesized in supernovae, C has to be formed mainly by other objects 
Unfortunately, because of their low metalicity, the past contribution of both C-rich 
gfants and WC stars to the elemental enrichment of these irregular galaxies has been 
small and thus these observations cannot be used to discriminate between these 
different types of stars. 


ra , te S' M ® r’ k P°' 2 corresponds to a C-enrichment 

timescale of the ISM of about 5x1 0 9 yrs. This may seem somewhat rapid, since the 
present C abundance in the ISM is very similar to the solar C abundance which 
presumably reflects the C abundance in the ISM 4.5x1 0 9 yrs ago However 
observations : allow an infal rate of unprocessed extragalactic material (eg., no heavy 
fh? wh.ch .sequal to the H I injection rate by stars (23) and this could increase 
the C-enrichment lifetime of the ISM by an order of magnitude. 


The last column in table 2 indicates the type of dust formed in the outflow of these 
SKf?!' (, he . best evidence for carbon stardust formation is for C-rich giants. Since 
about half of the carbon is tied up in gaseous CO in these objects, the carbon stardust 
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iniortinn rata is about 2x1 0* 6 yr 1 kpc' 2 . WC 8-10 stars also show evidence for 
drcurnsteha? cafbon sfardust. lire C fraction condensed in dust is however small 
?-3x1^2 461 The contribution of these objects to the carbon stardust budget of the 
n'a^axv (5x1 O' 8 MyrW 2 ) is thus minor. There is no conclusive eyidence for dust 
?orma\ion in supernova ejecta (47). Assuming that half of the carbon in supernovae 
rnSnses out yields 1 0’ 8 M n y r kpc' 2 for these objects. There is observational . 

223SS th S suqgests that s®ome novae form SiC as well as amorphous carbon grams 
AkhSi siHcates are not completely ruled out for the spectroscopic signatures 
observed we note that studies of meteoritic carbon also strondy suggest that novae 
?orm carbon stardust (cf., § VIII). We will assume that 25% of the C condenses in the 
n m d C staXst in novkl However, novae do not contribute much to the carbon 
d Sl budaet of the aalaxy. The total injection rate of C in solid form is then about 
SlO^M?yr^ kpc*, which is about 20% of the total carbon injection rate. This 
rather small fraction results from the large C contribution by objects which do not 
nTnduc? Carbon stardust (eg., O-rich giants, OB stars), the small fraction of C 
KnHonsPri^n the outflows from WC stars, and from the substantial fraction of C 
locked up in CO in C-rich giants. It should also be emphasized I that Ijhe dust 
supernovae and WC stars is ejected at high velocities (*1000 km s '). Appreciable 
desfruc°ion o1 this newly formed dust may take place during deceleration and mixing 
into the ISM due to sputtering and grain-grain collisions. 


C) The 13 C budget 

It is of some interest to examine the 13 C budget of the galaxy ^^^Here the situation 
iQ much less confused since He burning does not produce much C. As a result 
supernova^ and WC stars contribute negligibly to the Y3 C budget of the galaxy Again 
C-nch red giants can be very important. The 2 C/ 13 C ratio in those objects with the 

.IRC 1 is typically 45 

2^ translates then into a '^C mass-loss rate of 10 kpc Y r . • ^ r 
aiants and OB stars a solar ^ 2 C/' 3 C ratio is appropriate (=89j, which leads to an 
Fn^c\fonrete of3x10' 8 and 10* M 0 kpc' 2 yr^respectively. The ejecta from novae 
are exDected to be characterized by profoundly non-solar isotopic abundance patterns 
fofthe^iaht elements C, N and O (30,51). Novae may indeed be an important source of 
W ^ CaPculated r^C/^C ratios range frwn 0.1 to 0.5 (52,53). Adopting. 0.2 eads to a 
13 C mass-loss rate from novae of 2x10^ M @ kpc' yr . The total C injection rate 
is then about 3x1 O' 7 M 0 kpc' 2 yr 1 . 

This corresponds to an average 13 C mass abundance.in the ejecta of 1 .3x10**, which 
a cniar Annin infall of metal-poor extragalactic material seems to be 



is 

about 0.02 'which is about twice solar. As a result tne -u/ -n, ratio in mo .o.v. will 
slowlv (t*2x10'° yrs) increase compared to solar abundances. This comparison is 
independent of the (uncertain; infail rate of metal : poor extragalactic material. 

ActuaHy observations of the y2 C/“ 3 C ratio in the interstellar medium indicate an 
enrichment by about a factor 2 since the solar system formed (*4.5x1 0 3 yrs, 54) 
However since this enrichment timescale result from a comparison of two uncertain 
injection rates this difference with observations may not be too worrysome. 
observational studies of the 12 C/^ 3 C ratio as a function of galactic radius may be a 
very valuable tool in unraveling the contributions of different types of objects to the 
carbon budget of the galaxy. Preliminary observations indicate that there is no 
evident gradient in the ™CP 3 C ratio in tne solar neighborhood (55). This is somewhat 
Surprising in^ Sew of the large contribution of "pure^ 2 C by WC stars, which show a 
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marked J ncrease in surface density in the direction of the galactic center (42) This 
should be contrasted with the more even distribution of white dwarfs (eq red qiants 
descendants) and novae, which contribute most of the ^C. This might indicate that 

(wS»nHnl y = d nT nat ® ,h ,® buc !?®‘i 5 5>- ^ d0 note, however, that the 

O/H gradient in the galaxy is also fairly small (57) despite the large contribution 

from massive stars to the O budget and this argument may not be conclusive. 


V The Dust Condensation Zone 


We will consider three different types of objects, characterized as carbon-rich (eg., 

. 1 1 ). los |n g mass and showing evidence for dust formation. These are C-rich red 

» fnPiI2[SS B ° raalls S RCr B s ! a L s) ’ and 'ate-type, carbon-rich WoK Rayet sters 
Ww Wal J ftudied prototypes of these stellar classes are IRC 1 021 6 R Cr B 

of/n h.SUn Y? ?' 45, ^pectively. Each of these classes represents late stages in the 
f l X ol £ t J9 n ° f of different mass (and evolution). C-rich giants are low mass stars 
!mJhi^/nth^? h *l ave . be ®o (slightly) enriched in carbon through the dredge up of 
!y synthesfzed material. The evolutionary status of R Cr B stars is not well 
pnvpinio°H ' f jesumably, these low mass stars (=1 M 0 ) have lost most of their H 
thP^P^p?nnp onH P rocess,n 9 through the CNO cycle and possibly explosive ejection of 
fi®' r ? re n0W underway fr ° m a red supergiant stage to a white dwarf 

?t a 9®- T[j ese °A bje i ts are vary rare and have therefore been neglected in the C budget 
2 ,1. C V gala ? y - As dlsc ossed in 6 IV, Wolf-Rayet stars are an advanced staqe in the 

fhe surfarp Tahip'?!?? ? a - S \ M>25 )■ where He burnin g products are exposed at 

tne surface. Table 3 lists typical parameters for these prototypes and their H 

58 C 63) S |t e ^h22iH V hp£ e ' T h hese hav ® been C0 . m Pile d from the recent literature (41-43, 
d be . em P basized that most of the parameters reported actually show a 

Srther' dIriftS W " hm 6a0h ClaSS ' A few °< the reported Parameters neeS tome 
A) The photosphere 

Thenar radii have been calculated from the luminosities assuming that the stellar 
n^/dl^iT/pT 1 b 2 e re P reser) ted by a black body at the effective temperature, T pff ?R? = 

[L /4jca] T e « 2 , where a is Stefan-Boltzman’s constant). The total pressure ?n the 

t P e °Q cnante ’it has [Sn S tak th0 f deplh W wP re the continuum emission originates. For 
Qtprc iFhie S baen t taken fron ? model calculations (65) for log g=1 . For the R Cr B 
atara '- has been ® st| mated assuming hydrostatic equilibrium and the Rosseland mean 
opacities appropriate for these starsM The total pressure in the photosphere \T 
meaningful concept for Wolf-Rayet stars where the continuum emission 

th is wavTachSl^f^Ifh well i n j£ the °- utflo V v and ’ thus > the stel,ar radius (defined 
tn^ way) actually falls beyond the sonic point of the outflow 


hm?^!?v?t'2f ed *9 / ^2 d C L ° ratios for an optically selected sample of C-rich qiants are 
mLTloss^a^ whirh a ha 1 " k’ res P e ctively (67). However, for our purposesttie high 
Sto S riJ? S Th S i arS ’ T2i ch , have been selec ted against in this sample, are of more y 
ah2SfinI^ e 3 uate ^ y aluas seem reasonable and are also in agreement with 
abundances derived for planetary nebulae (58,67). The ^C/I3g ratio d eDenc | s ctronnlv 

on the amounts of He burning shell and equilibrium CNO cycle material mixed into th 9 e V 
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Table 3: Physical Conditions in the Carbon Dust Condensation Zone 


Type 


C-rich Giant R CrB 


WC 8-10 


STELLAR PARAMETERS 


- 2 ) 


L* (L J 

Tel/) 

R*(RJ 
P t0 ,(dyn cm 
Period (days) 
C/H 
C/He 
C/O 

13 C/ 12 C 


10 4 

10 4 

10 5 

2500 

7000 

20000 

530 

67 

25 

4(3) 

3(3) 

” 

200 

40 

““ 

6 (-4) 

300 


7 (-3) 

3 (-2) 

8 (-2) 

1-2 . 

2 


3 (-2) b 

<2 (-2) 

|0] a 


CIRCUMSTELLAR SHELL: CONDITIONS AT THE CONDENSATION RADIUS 


M 


.-1 


,of e yf 
(km s 1 ) 

R c (R.) 

T d (K) 

T g (K) 

P, ot (dyn cm 


- 2 \ 


carbon 


(” ”) 


1 exi 


p/ T col 


observed 

composition 


1 (-5) 

15 

5 

=1000 

=800 

6 (-5) 

4 (-8) 

500 

molecular 
(see table 4) 


1 (- 6 ) 
200 
10 

=1000 


3 (-6)® 

1 <; 7 ) 

40° 

c, c 2 ,cn 


3 (-5) 
1500 
500 
=1000 


3 (-8) c 
2 (-9) c 
2 C 

C + . C 2+ , 

r>3+ r 


notes: 


a) These abundance ratios refer to theoretical calculations rather than observations. 

b) This ratio varies widely from object to object within the range 3.5-1 00. 

c) These numbers have been calculated assuming a spherically symmetric outflow and the actual 

densities and pressures may be much higher (see text). 


atmosDhere bv the dredge up process during the helium shell flash at the tip of the 
asymptotic giant branch(68). Observed values range from 3.5 to 100 (ie., .solar). The 
value shown for IRC 10216 (cf., table 3) seems to be typical for high mass-loss 
C-rich giants (49). 
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THE STRUCTURE OF C-RICH RED GIANTS 


* 


3 AU 


♦ 


♦ 


10 4 AU 


♦ 





Figure 4- a schematic representation of the structure of red giants (top) and their surrounding 
circumstellar shells (bottom). 1 AU is the sun-earth distanced .5x1 0* 3 cm. For clarity the various 
inner zones are not drawn to scale. A star on the asymptotic giant branch consists of a central, 
carbon-oxygen, white dwarf (=1 0 3 AU) surrounded by a thin shell - in which nuclear burning takes 
place - and an extensive, convective envelope. Strong shock waves in the photosphere transport matter to 
large distances (=5R») where dust condensation takes place (stationary layer). Radiation pressure on the 
dust accelerates the dust outwards and the gas is dragged along. 



B) Structure of the outflow 

A star on the asymptotic giant branch (AGB) has exhausted its hydrogen and helium 
nucleosynthetic fuel in the core, but is not massive enough (<8 M @ ) to explode as a 
supernova (type II). The core then collapses into a carbon-oxygen white dwarf (=10 
cm) which is surrounded by a nuclear active zone in which alternately hydrogen and 
helium burns (cf., fig. 4). This core and nuclear active zone are embedded in an 
extensive (=5x1 0 1 3 cm) convective envelope (68). The envelope is unstable and 
Dulsations drive strong shock waves into the photosphere, transporting matter to 
arge distances from the stellar surface. This material forms a (quasi) stationary 
ayer at about 5R* (=2x1 0** cm), where solid dust particles condense put. Radiation 
pressure on these dust grains leads then to the formation of an extensive (=10 cm), 
expanding (=10 km s' 1 ) dust shell. The gas is collisionally coupled to the dust and is 
dragged along. During the last part of this evolutionary phase ofa 'jw- mass star 
when the mass of the convective envelope becomes small (=1 0 M Q ) a fast, low 
density wind (ie., "superwind") from the white dwarf compresses the surrounding gas 
and dust shell, which is ionized by UV photons from the white dwarf. The AGB star and 
its surrounding shell has become a planetary nebula. The nuclear fuel is exhausted and 
the white dwarf will eventually cool into oblivion. During the penultimate phase, 
some Dlanetary nebula nuclei develope atmospheric characteristics of carbon-rich 
Wolf-Rayet stars (ie, population II WC stars; 69). 

The mass-loss of population I Wolf-Rayet stars (ie., M*>20 M @ ) is due to radiation 
pressure on atomic lines. Dust forms only at large distances frorn the stellar surface 
and in contrast to red giants, plays no dynamic role in the outflow. These radiation 
driven winds from Wolf-Rayet stars are highly unstable and strong shocks are 
expected to develop (42). Velocity amplitudes may be as large as 500-1 000 km s . 

For dust formation, the most important effect of these instabilities is perhaps the 
development of dense pockets of gas which can shield the budding molecule and soot 
formation processes from the disruptive effects of stellar and 

photons. Indeed, the small fraction of C condensed in the form of dust and the absence 
of IR variability suggest that dust formation occurs in a large number of highly 
localized dense regions in the outflow. The increase in dust mass towards later-type 
WC star is then attributed to a larger number of such regions, caused perhaps by a 
larger susceptibility to this instability. 

In R Cr B stars, ejection from the pulsating stellar photosphere seems to take place 
in periodic, extremely dense, but localized, "puffs". During the expansion the 
temperature in these ejected gas blobs drops and dust is formed (62). Due to the 
localized nature of the ejection process, the stellar disk is only sporadically covered, 
aivinq rise to ocassional deep minima in the observed light cun/es of these objects. 
Dust condensation takes therefore place at much higher densities than derived from 
the nominal mass-loss rates assuming spherical outflow at the terminal velocity. 

C) Physical conditions in the dust condensation zone 

IR observations of the vibrational transitions of CO show that C-dust condensation in 
red qiants takes place at a relatively large distance from the stellar photosphere 
(=5R*- 63) This is in agreement with IR speckle measurements, which show similar 
mnner shell radii (58). For R Cr B and WC stars the condensation radius has been 
estimated from the condensation temperature assuming reasonable absorption 
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properties for the dust (64,70). The derived inner radius for the WC stars is in good 
agreement with IR speckle studies of Ve 2-45. For red giants the gas temperature at 
the condensation radius is measured to be about 800 K (63). No measurements for WC 
and R Cr B stars are available. The total and partial carbon pressures, P tnt and P r at 
the cpndensation radius have been determined from the mass-loss rate, tne outflow 
velocity (eg., P^tcR^v^]) and the (gas or dust) temperature. As remarked above, 
because of clumping, these values may not be very relevant for R Cr B and WC stars 


It is obvious from table 3 that carbon dust formation occurs under a wide variety of 
conditions. The effective temperature of the central star, for example, ranges from 
2500 to 20000 K and this will have a profound influence on the chemical composition 
of the expanding shell (cf., below). The elemental composition of the outflow also 
varies dramatically from mainly hydrogen to mainly helium and carbon in the three 
classes of objects. However, dust formation always seems to occur at a condensation 
temperature of about 1 000 K. An important parameter to consider in dust formation 
studies is the ratio of the expansion timescale, t 6 „ d> to the collision timescale, T rn , 

(71 ). The expansion timescale is approximately given by R Jv^ (cf., table 3). 
Anticipating the discussion in § VI, dust formation is initiated by collisions with the 
most abundant, carbon bearing, gas phase molecule. For C-rich giants this is acetylene 
and, assuming a typical neutral-radical reaction rate (2x1 0' 11 cm* s' 1 ) the 
calculated ratio is embarrassingly small (cf., table 3). This implies that kinetic 
effects (eg., activation barriers) are extremely important and that the dust formation 
process may be far from thermodynamic equilibrium. Essentially, the dust formation 
process has to be extremely efficient and shpuld occur upon collision. Although grain 
growth probably occurs at this rate, such a high efficiency is difficult to understand 
for the initial nucleation step, which is probably inhibited by considerable activation 
barriers. On the other hand formation of the initial condensation nuclei (eg., PAHs- see 
§ VI) in or slightly above the stellar photosphere, will be inhibited by the strong 
shock waves traversing this region (72). This presents an enigma for any stardust 
formation theory. It has been argued that the observation of a stationary layer and the 
IR speckle results are perhaps not in direct conflict with dust nucleation starting in 
the photosphere (58). Indeed, small PAH nuclei possess a large electronic band gap 
(cf., § III) and will not experience much radiation pressure until they have grown to 
about 50 A. Their dynamic influence on the outflow is therefore small. Likewise, if 
only a small fraction of the carbon is locked up in these species, they will not be 
detected in the IR (unless a bright FUV source is present as in planetary nebulae) The 
observations do imply, however, that further growth to dust particles is inhibited 
until about 5R* above the photosphere. Since nucleation is the bottleneck to dust 
formation (further growth of nuclei into large dust grains is expected to be rapid) 
these observations do indirectly constrain nucleation to the stationary layer. 

D) Chemical composition 

The last row of table 3 summarizes the chemical composition in the outflows from 
these objects. For R Cr B stars only a few very simple carbon-bearing species have 
been observed. The winds from WC stars are completely ionized. The abundance of 
neutral carbon is about 1 0' b of the elemental carbon abundance in these stars A wide 
variety of molecular species have been observed in C-rich late type giants. These are 
summarized in table 4. This data again refers specifically to the outflow from the 
well-studied C-rich giant, IRC 10216 (73-76). The IR observations refer to pencil 
beam observations against the IR stellar and dust continuum and are thus heavily 
weighted towards the high density region close to the star. Molecular abundances can 
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Table 4: Molecular Composition of Outflow from C-rich Giants 


I. Infrared Observations 3 

Species Nj 

[cnrn 


II. Radio Observations* 5 

Species Nj 

[cm 2 ] 


CO 

2.0(20) 

6(-4) 

HCN 

1.(16) 

CpHp 

HON 

3.0(1 8) c 
>1.5(18). 

9(-6) 

>5(-6) 

HCoN 

hc 5 n 

3.(15) 

6.(14) 

NHo 

1.0(17)° 

3(-7) 

hc 7 n 

3.(14) 

CH 4 

1.8(16) 

5(-8) 

HCgN 

2.(14) 

SiH 4 

C|H 4 

4.0(15) 

4.0(15) 

6.0(19) 

2(-7) 

M-8) 

2(-4) 

HCi i N 
CHqCN 

cr? 

c 3 h 

2.(14) 

4.(12) 

9.(15) 

1.(14) 




c 4 h 

3.(15) 




c 5 h 

1.(14) 




CgH 

2.(14) 




c 3 n 

1.(15) 




c 3 h 2 

2.(14) 




sic 2 

2.(15) 




SiO 

1.(16) 




SiS 

7.(15) 




HpS 

1.(13) 




cs 

4.(15) 




C 2 S 

1.(14) 




C3S 

1.(14) 




NaCI 

5.(12) 




AICI 

5.(12)' 




KCI 

5.(12)' 


n o( es - 

a) Data taken from the compiliation in reference (73). Abundances calculated assuming a CO/H 2 

abundance of 6x1 0' . 

b) Data taken from references (74-76). 

c) 600K C 2 H 2 With N= 2.0(1 9) also present. 

d) Cooler (<200K) NH 3 also observed. 

e) Column densities assumed to be equal to that of NaCI. 


be derived from the measured column densities, assuming a size scale (ie., the inner 
radius). The radio studies, in contrast, are dominated by the cooler, less dense 
material at large distances from the star. Since the observed intensity measures the 
total emiting volume, the column densities derived from the radio observations 
depend critically on the assumed (ie., unknown) inner radius (74). No direct 
comparison is therefore possible between the two sets of column densities given in 
table 4. 
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OMPOSITION OF THE STATIONARY LAYER 



Figure 5- The chemical composition of the stationary layer calculated assuming thermodynamic 
equilibrium (adapted from reference 63 ). 
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A strikingly rich variety of molecules has been observed in IRC 10216, much more 
than in any other circumstellar shell. To a large extent, this is probably due to its 
proximity to the sun?d«200pc) and other C-rich envelopes are likely to be as 
chemically diverse. The molecules in table 4 can be divided in four categories: 1 ) 

Simple and stable molecules expected in a C-rich environment (eg., C 2 Hp> i HCN CH,, 

Co Ha - see below). 2) A series of acetylenic chains (cyanopolyynes; eg., hC 3 N, HC 5 n). 3) 
A series of reactive radicals (eg., C 3 H, C 4 H). 4) A set of molecules involving S, Si and 
Cl This diversity reflects the large number of chemical processes, which are of 
importance in circumstellar envelopes. These include (nearly) thermal equilibrium 
processes in the stellar photosphere and inner envelope, shock processes, UV 
photochemical processes, grain condensation and grain surface reactions (61). All of 
these can influence the chemical composition of the dust condensation zone. Figure 5 
shows the abundance of some simple carbon bearing molecules as a function of aqs 
temperature calculated assuming thermal equilibrium at a total gas pressure of 1 0 ' 3 
dyn cm^ and a C/O ratio of 2 (63). This pressure is perhaps a little higher than 
expected for the dust condensation zone in C-rich giants (cf., table 3). At high 
temperatures, most of the carbon is locked up in CO and C 2 H 2 - At lower temperatures, 
CH,, becomes important and even some of the CO is converted, which results also in an 
increase of the abundance of other O-bearing molecules (eg., SiO, H 2 0; 63). Besides 
acetylene attention should be drawn to the relatively high abundance of CoH at high 
temperatures, which may play a key role in carbon dust condensation (cf., s VI). Note 
that these calculations ignore the dust condensation process. Similar calculations 
have been performed including dust condensation (again assuming thermal 
equilibrium), but for a high pressure ( P tot = 10 3 dyn cm’ 2 ) more appropriate for the 
stellar photosphere rather than the condensation zone (77). Moreover all of these 
calculations ignore chemical kinetics which 

certainly plays an important role in determinig the molecular composition and 
ultimately carbon dust formation (cf., § VI). Finally, atomic H may also play an 
important role in soot formation. H 2 is the most abundant molecule in these 
equilibrium calculations and atomic H is a trace species. However, shocks driven into 
the extended photosphere by the stellar pulsations in red giants may actually lead to 
Ho dissociation and the H/H 2 ratio may be far from equilibrium since H 2 reformation 
wnl be slow under these conditions (78). 


VI Chemical Pathways to Carbon Stardust 


In this section we will discuss the formation of carbon dust in the outflow from 
carbon-rich objects. We will concentrate on C-rich giants, since they dominate the 
carbon stardust budget (cf., § IV) and since the physical conditions in their outflows 
are much better determined than those in WC stars or R Cr B stars. Moreover, their 
chemical composition resembles that of sooting flames and we can draw on an 
immense chemical literature on soot formation under such conditions (79). Some of 
the differences that can be expected to be important in other objects will be pointed 
out. 

A) Thermodynamics and carbon condensation 

Most studies of stardust formation are based on thermochemistiy, which considers 
the condensation of solids out of a (slowly) cooling gas that originally contained no 
solids (71,80). The "condensation" sequence is then determined based on 
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THERMODYNAMICS OF CARBON DUST FORMATION 


P, atm 

1(T 10 io~ 5 



Figure 6 . The phase diagram appropriate for C-rich red giants (adapted from reference 80). The line 
labeled solid carbon encloses the regime in which graphite is stable. Surrounding it, the position of the 
various, stable, carbon-bearing gas phase molecules are shown. To the left (right) of the dashed curve 
hydrogen is mainly atomic (molecular). The solid line in the lower right hand corner deliniates the 
transition from the CO to the H 2 0 dominated regime. The conditions appropriate for the photospheres and 
dust condensation zones of C-rich giants are also indicated. Due to the presence of strong shocks, the 
actual path of a cooling gas element in the outflow is far more complex than indicated. See text for details. 
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thermodynamic equilibrium. Essentially, the free energy. of many atoms, molecules 
and solids are used to derive their relative abundances in equilibrium. Condensation 
sequences have been published for many different pressure, temperature, and 
elemental abundance regimes (71 ,80-82). Figure 6 shows the result for a C-rich (and 
H-rich) atmosphere (80). The thick solid curve encloses the regime where graphite is 
stable. Outside of it gas phase carbon species are more stable. Their relative 
abundance depends on the total (H) pressure. At the lowest pressures C, C 2 and C 3 are 
most important Between 10' 9 and 10 4 dyn cm' 2 acetylene dominates, while at high 
pressures CH, ties up the carbon (80). To the left (right) of the dotted curve atomic H 
(Ho) is most abundant. The conditions appropriate for the photosphere and the 
condensation zone in C-rich giants are also schematically indicated (cf., table 3). As 
expected from observations (59) the conditions in the stellar photosphere correspond 
to the acetylene dominated regime. Carbon stardust formation occurs in a regime 
where graphite is thermodynamically favored (cf., fig. 6). The long-dashed line 
schematically indicates the path of a cooling and expanding gas element. The actual 
oath is however, much more complex, due to the strong shocks traversing the 
extended photosphere (cf., § VB). Nevertheless, condensation does not seem to take 
place immediately upon passing the solid carbon stability line. The flow is, thus 
hiahlv supersaturated when the grains actually condense and an amorphous rather 
than crystalline material will form (cf., 1). The thermodynamic stability curve for 
such a material is, however, expected to be very similar to that of graphite (83). This 
late onset of condensation may reflect the importance of kinetics (eg., 
time-dependent nucleation) as well as other non-equilibrium effects. 


The actual grain growth process is generally calculated using classical homogeneous 
nucleation theory (80,84). This assumes that the abundance of a cluster containing N 
atoms can be derived from the monomer concentration using the free energy of this 
cluster The latter is often derived from the experimentally measured, bulk surface 
tension (84), although this may lead to large errors for small clusters (85). For small 
clusters, the surface free energy makes alarge contributions to the total free energy. 
As a result, below a critical cluster size (which depends on the degree of 
supersaturation), clusters are more likely to evaporate than to grow. Statistical 
fluctuations are then responsible for producing clusters larger than this critical size, 
which can then grow on to actual dust grains (83,84). The formation of the critical 
cluster forms thus a bottleneck in the dust gram nuc eation. Typical the final number 
of monomers in a grain is then approximately given by h/lnfp)] 1 '' 3 (80), where p is 
the number of monomer collisions within a relevant timescale (ie.* T ex yT C q|). hor the 
conditions in C-rich red giants (cf., table 3) this resu ts in N«5x10 5 or ^ particle size 
of about 1 0OA. Since the initial steps are expected to be associated with considerable 
activation barriers, the actual grain size will be much smaller. This sma I expected 
grain size reemphasizes the dust nucleation problem alluded to in the previous 
section. 


Chemical reaction rates are, to a large extent, governed by kinetic rather than 
thermodynamic factors, a point whicn has largely been ignored in stardust 
condensation studies. In particular, nucleation theory was developed for the 
formation and growth of liquid droplets where the bulk surface tension is more 
accurately known than for solid clusters. Indeed, the concept of a surface (or 
surface tension) is not even well determined for small molecular clusters (85) and 
such reasoning can easily lead to erronous results. For example, it has been argued 
that, if the surface free energy of tetrahedrally-bonded carbon clusters is equal to 
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figure 7- Transmission Electron Microscope photographs of a typical amorphous carbon sample produced by 
striking an arc between two carbon electrodes (courtesy of reference 92). Top: Carbon necklace. Bottom: a 
single grain. 
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that of aromatically-bonded clusters, diamond clusters would be favored over 
graphitic clusters in the carbon condensation route around stars, essentially due to 
the difference in volume density (86). However, aromatic clusters grow at the 
periphery of the aromatic planes and the surface free energy is associated with the 
dangling bonds of these peripheral carbon atoms. In contrast, diamond clusters grow 
three dimensionally. As a result the number of "surface" carbon atoms and their 
associated dangling bonds per carbon atom in the cluster is much larger in a 
"diamond" cluster than in a '’graphite" cluster. It is only at high pressures, where 
volume plays a dominant role, that the high volume density of diamond favors 
diamonds over graphite. As a result aromatic clusters (eg., PAHs, see below) are 
expected to form around stars in agreement with astronomical observations and with 
laboratory studies of soot formation. 


Althouqh one may attempt to circumvent this problem in the definition of the cluster 
surface free energy in one way or another (84), the application of classical nucieation 
theory to stardust condensation is even more fundamentally flawed. In particular the 
monomers (09., atomic C) out of which graphite or amorphous carbon has to be formed 
do not even exist, as such, in the outflow from C-rich giants where most of the 
elemental carbon is in the form of CO and C 0 H 2 (87). A similar argument holds also 
for the condensation of silicates around O-rich giants. For a meaningful assesment of 
carbon dust formation, the chemical pathways which convert these molecules into 
solid dust grains has to be identified (79,88). 

B) Soot formation 

The processes that convert C 2 H 2 into carbon grains in the outflow from C-rich red 
giants are probably very similar fb those occuring during the gasphase pyrolysis ot 
hydrocarbon molecules. There is an extensive literature on soot formation in 
combustion environments, in a large part driven by concerns on internal combustion 
engine efficiency, pollution, and healtn hazards (ct., 89-91 and references therein). 
Here we will briefly summarize the salient points. 


Soot is a general side product of the combustion and pyrolysis of hydrocarbons. 
Morphologically, soot is independent of generation method and fuel type. Transmission 
electron microscopy shows that soot particles consist of a large number of small 
(1 00-1 OOOA) spherical partic es arranged like pearls on a chain. Figure 7 shows a 
typical example generated by striking an arc between two amorphous carbon 
electrodes (taken from 92). Similar results are obtained when burning hydrocarbons 
such as methane, acetylene or benzene in air at room temperature and in pyrolysis 
studies in shock tubes and flow reactors (89-95). These particles are made up of 
microcrystallites with sizes in the range 10 to 25A. X-ray analysis shows that the 
carbon atoms are arranged in a graphitic planar structure and thus such 
microcrystal I ites resemble large PAH molecules. Indeed, smaller members of the 
condensed aromatic ring compounds (eg., PAHs) are readily extracted from soot 
particles using organic solvents. Although the microcrystalline planes are parallel to 
each other, they are more randomly stacked in soot particles than in graphite (95). 
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Figure 8 . The structure of the various molecules and particles involved in one possible carbon stardust 
formation scheme. Acetylene is the precursor gas phase molecule. Various chemical processes convert 
these molecules into benzene molecules. Further chemical growth involving acetylene leads to the 
formation of large polycyclic aromatic hydrocarbons. Several of these planar molecules stacked together 
form a platelet. These platelets are the building blocks of soot particles. These platelets as well as the 
layers within them are generally cross-linked by tetrahedrally bonded carbon atoms and chains (see text 
for details). 
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The chemical nature of the soot can vary within a flame. When first formed, the 
particles contain a large number of hydrogens (C/H=1 , when a=50A), which are 
distributed in between the aromatic layers. ESn studies show that these small 
particles have a pronounced free radical character and they continue to grow through 
surface reactions. Both the hydrogen content (C/H-0.1) and radical nature decreases 
during this growth. At a size of about 200 A rapid surface growth stops and the 
particles aggregate in larger grains. 


Three distinct steps can thus be discerned in the formation of soot (89-91): 1) the 
breakdown of fuel molecules to precursor molecules which react rapidly to form 
larger species; 2) chemical growth that results in the formation of numerous small 
particles (a>1 0A); and 3) particle growth by surface reactions as well as by 
clustering and agglomeration. In flames these processes compete with oxidation 
processes. In the outflow from C-rich red giants, however, the C/O ratio is larger 
than 1 and essentially all of the oxygen is tied up in CO. In other words, all of the 
available carbon has already been burned in the stellar photosphere and the conditions 
resemble those of a fuel-rich flame. Thus, while the incandescent soot particles in 
flames, which give rise to its characteristic yellow color, will typically burn up 
quickly in the air at the top of the flame, this oxidation process will be of no 
importance in C-rich giants. 


C) PAHs and the formation of carbon stardust 

The critical intermediates in the carbon condensation route have to be highly stable 
against dissociation and yet posses a high reactivity towards polymerization. Several 
different chemical pathways to soot formation have been proposed based on neutral 
radicals, ions, PAHs, polacetylenic chains, or Buckminsterfullerene as intermediaries 
(89,90,96,97). The high tendency of aromatic fuels to soot, the high stability of PAHs 
as well as their structural similarity (the aromatic C backbone) has led to the notion 
that PAHs are the building blocks of soot particles. In the case of C-stardust 
formation around C-rich giants, the importance of PAHs is underscored by their 
ubiquitous presence in C-rich planetary nebula (7). Following references (79,88), we 
will identify PAHs as the key intermediates in C-stardust formation in these objects. 

The basic steps in soot formation in this scheme are schematically shown in figure 8. 

As remarked above, other pathways have been proposed in which PAHs are unreactive 
sidebranches (ie., dead-ends), and the following discussion should be considered 
illustrative rather than authorative. Moreover, as will be discussed in 6 E, it is likely 
that in different environments different mechanisms dominate C-stardust formation. 


In the outflow from carbon-rich giants, acetylene and its radical derivatives are 
likely to be the precursor molecules from which soot is formed and the soot-limiting 
step is the formation of the first aromatic ring (cf., fig. 8). Rapid chemical growth of 
this ring then forms larger PAHs. At some stage clustering of PAHs becomes 
important and small platelets are formed in which small PAHs are randomly stacked. 
Simultaneously, further chemical growth will take place, in particular at the 
periphery of the aromatic planes. This may lead to cross-linking of the planes by 
tetrahedrally bonded carbon. Some crosslinking may also occur randomly dispersed in 
the aromatic planes. Finally, these platelets can cluster to form a spherical soot 
particle (cf., fig 8). In flames, these spherical soot particles cluster further to form a 
^pearl-like" necklace (cf., fig. 7). 
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THE CHEMICAL PATHWAY FOR CARBON DUST FORMATION 



fi gure 9 - 0r) e possible chemical pathway towards soot (adapted from reference 103). Top: the formation of 
the first aromatic ring, starting from acetylene. Bottom: further chemical growth takes place through 
alternating steps of radical site formation and acetylene addition. Branching reactions have been omitted 
for clarity (see text for details). 
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The temperature and pressure history of a parcel of gas in the outflow of a 
carbon-rich giant is very complex due to the strong shock waves propagating through 
the extended photosphere which lead to large temperature and density excursions. 
Dust formation will then occur far above the photosphere (R=5R*) where the activity 
generated by the shock waves has died down (72). PAH formation in flames starts 
typically when the temperature drops below about 2000K (98). The temperature of the 
extended red-giant photosphere is not determined by radiative equilibrium but by 
shock passage. A postshock temperature of 2000K corresponds to a shock velocity of 
a few km s* ' and shock calculations in Mira atmospheres indeed indicate a typical 
distance of about 5R* above the photosphere (99-102). 


Figure 9 shows the first steps in one possible chemical pathway to soot derived from 
detailed kinetic model studies of shock-tube pyrolysis of acetylene at about 1700 K 
(78, 79, 103). A very similar cyclization model has been proposed to explain measured 
product distributions in sooting laminar flat flames with a variety of hydrocarbons 
(104). Many different branching reaction pathways are possible at all of the 
intermediate steps. For clarity, these have been ignored, since theoretical studies 
indicate that they are of lesser importance. There are two steps in this reaction 
scheme: a series of reactions describing the formation of the first ring out of 
acetylene and one describing further chemical growth. Formation of the first 
aromatic ring forms a bottleneck for the soot yield. It is started through the 
formation of the ethynyl radical (C 2 H) through H abstraction. This radical can then 
attack the triple bond in another acetylene molecule and react to form C 4 H 5 . Two 
steps of H abstraction by H atoms, followed by reaction with yet another acetylene 
yields C 6 H 5 . Cyclization of this radical via interaction of the unpaired electron with 
the triple bond leads then to the first aromatic ring. Once the first ring is formed 
subsequent chemical growth consists of (more or less) alternating steps of the 
formation of a radical site through H atom abstraction, acetylene addition and 
cyclization. This is illustrated for the formation of the small PAH, naphthalene. 

Further chemical growth will take place through the fused PAHs pyrene, coronene, and 
ovalene (105). 


The stable molecular intermediates, including phenylacetylene (C 4 H 4 ), acenaphthalene 
(C 12 H 8 ) and the pericondensed PAH pyrene (C 16 H 10 ), present in this reaction scheme 
play an essential role. The reactions leading from one radical to another are 
reversible at these temperatures and quickly set up an equilibrium distribution. The 
reactions leading to these stable intermediates, however, are irreversible and thus 
"pull" the reaction sequence towards molecular complexity and soot formation. The 
importance of pericondensed aromatic hydrocarbons, such as pyrene, coronene and 
ovalene, in this scheme rest mainly in their high thermodynamic stability, resulting 
from their maximum resonance energy and minimal stress energy (105). The enthalpy 
of formation can also be maximized by minimizing the entropy change by returning as 
many H 2 molecules as possible to the gas phase. This favors molecules with an as 
short as possible boundary structure and may actually lead to (smaller) aromatic 
molecules having pentagons or squares in their carbon skeleton (78). The molecule 
cyclopentanaphthalene is an example of this latter principle. The stresses due to the 
presence of pentagons in the carbon skeleton will lead to warping of the intrinsically 
planar structure of PAHs. If present in an ordered fashion this can even result in 
spheroidal molecules such as Buckminsterfullerene (C 60 ; 106), which are closed upon 
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Table 5: The most abundant PAHs in flames 9 


Species 


formula 


abundance* 5 


Naphthalene 

ClO H 8 

100 

Acenaphthalene 

Ci 2 h 8 

100 

Phenanthrene 

C14H10 

30 

Pyrene 

Cl 6 H 10 

30 

Benzo[g,h,i]fluoroanthene 

^18^10 

10 

1,12-Benzperylene 

C 22 Hi2 

3 

Coronene 

Soot 0 

C 24 H 12 

1000 d 


notes: 

a) Taken from (98) and references therein. 

b) Mass abundances in the burned gas normalized to coronene. 

c) Typical soot diameter =250A, corresponding to about 10 5 C atoms. 

d) Measured when chemical growth has stopped. However, further physical growth (ie., clustering) 
may increase this ratio further. 

themselves and do not contain any hydrogen. When present in a disordered fashion, 
pentagons will lead to a curling of the aromatic planes without actual reconnection 
(96,9/). However, cyclopenta groups are relatively unstable and may be partially lost 
along the reaction pathway (103). Whether they will persist in the carbon 
condensation route probably depends on the relative importance of thermodynamic and 
kinetic considerations. 


Table 5 list some typical PAHs detected in sooting flames with an estimate of their 
relative abundance (98). In the case of well-controlled combustion mainly the 
thermodynamically most stable PAHs are formed and survive in the burned gas. 
Although the absolute amount of PAHs formed depends critically on the physical and 
chemical conditions, the relative concentrations of these stable PAHs do not. This 
table is mainly concerned with the pericondensed PAHs which are thermodynamically 
favored (105). The abundance of hydrocarbons with the same number of carbon atoms 
but a larger hydrogen content and a less stable structure is much lower (98). For 
example, the abundance of the 5 ring, benzo[g,h,ilfluoroanthene and 
cyclopenta[c,d]pyrene is two orders of magnitude larger than that of triphenylene, 
chrysene, and benzoanthracene. Likewise, apart from coronene (C 24 H 10 ), no isomers of 
this fully condensed benzonoid have been detected in sooting flames. Finally, the high 
abundance of species containing pentagons (eg., acenaphthalene) is not surprising in 
light of the discussion above. Searches for these PAHs in the outflow from C-rich 
grants will be very valuable to test the similarity of dust formation around C-rich 
giants with that of sooting fjames, to elucidate the actual carbon condensation 
pathway and to determine the relative importance of different interstellar PAH 
species. Of course, the distribution of PAHs ejected by these objects will be further 
modified by processes in the diffuse interstellar medium. In particular, UV 
photodissociation and photoisomerization will play an important role in weeding out 
less stable as well as small PAHs in the diffuse interstellar medium (107). 
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The formation of carbon stardust in the outflow from carbon-rich red giants is, thus, 
expected to be alike to the formation of soot in flames. Numerous laboratory 
experiments suggest that PAHs are, if not an essential intermediate, at least an 
abundant side product of the soot formation process. The infrared spectra observed 
towards many planetary nebulae (the red giants "daughter") show emission features at 
3.3, 6.2, 7.7, 1 1 .3 pm as well as broad underlying emission component (7). Although 
the’observed spectra compare very well to that of carbon soot particles, energetic 
considerations imply that tne species emiting the emission features are very small 
(20-50 carbon atoms; 10,1 1). they are, thus, a collection of free-flying molecules 
rather than clustered in soot particles. The broader emission components as well as 
the far infrared continuum is evidence for the presence of larger PAHs and soot 
particles (58 108). Finally, we note that the ratio of C in PAHs relative to that in 
soot observed in planetary nebulae (=20; 108) is in reasonable agreement with that 
measured for sooting flames (cf., table 5). 

D) Non-LTE and the formation of C stardust 

The conditions in the outflow from carbon-rich objects are probably far from being in 
thermodynamic equilibrium and this should have a profound influence on the carbon 
dust condensation process. For example, the vibrational excitation temperature of 
molecules in the dust condensation zone may be higher than the kinetic gas 
temperature due to radiative excitation by stellar photons or the trapping of line 
radiation (109). Since reaction rate constants are often very sensitive to vibrational 
excitation, some reactions may proceed much faster than expected from the kinetic 
gas temperature. The chemical composition may also be far from equilibrium. This is 
m particular true for the expected H/H 2 ratio. In the photosphere H 2 is formed by 
three body reactions and tne H/H 2 ratio freezes out when the density drops below 
about 10' ' cm' 3 (110). This results in a disproportionally large abundance of atomic 
H farther out in the flow. The atomic H abundance is also influenced by the strong 
shock waves in the extended photosphere. Shock velocities can be as high as 25 km 
s' 1 and Ho dissociation will be important. Since the post shock density may be less 
than the freeze-out density, large deviations from the expected H/H 2 ratio may occur. 


Finally, ionization may play an important role. Although no ionizing radiation is 
expected from a cool red giant, chromospheric radiation may play a (ionizing) role in 
supergiants. Moreover, chemiionization has been shown to be an important ionizing 
process in sooting flames. Of particular importance are the reactions 

O + CH ~> HCO + + e‘ 

CH(A 2 A or B 2 2) + C 2 H 2 --> C 3 H 3 + + e‘ 

These ions are present in the reaction zone in sooting flames and their abundance 
falls off in the burnt gasses. Reactions between these ions and neutral molecules are 
rapid, as is isomerization into their most stable form (eg., PAH ions for the larger 
ones). They may actually play an important role in the soot nucleation and 
condensation process (90, 1 1 1). As for the neutral PAHs, pericondensed structures, 
such as C 13 Hq + , C^yH^ 4 " and C 19 H 11 + , are thermodynamically favored and observed to 
be most abundant in flames. In the circumstellar case, however, oxidation occurs in 
the stellar photosphere and little atomic oxygen (or CH) is expected to be present in 
the outflow. Cosmic rays which have sometimes been considered as an ionizing agent 
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are unlikely to play a role, since they are tied to the galactic field (gyro radius = 

1 0' J cm), which will be wrapped around the bubble blown by the stellar wind. Thus, 
ionization processes may be much less important in circumstellar shells around 
carbon-rich giants than in sooting flames. 

E) Carbon stardust formation around WC 8-10 and R Cr B stars 

The formation of carbon dust around R Cr B stars and WC 8-10 stars is expected to 
follow a completely different chemical pathway than in the outflow from C-rich 
giants. These objects have lost most of their H-rich envelope either through nuclear 
burning or through previous mass-loss. The reaction scheme derived above from 
experimental studies of the pyrolysis of hydrocarbons in sooting flames (cf., figure 
8), which may be very relevant in the outflow from C-rich red giants, will play no 
role in the carbon dust formation in these objects, since no C 2 H 2 , C 2 H or H is present 
Carbon soot formation is now expected to be initiated through the formation of 
polyynes (eg., acetylene-like C chains withouth H; ie., C 2 , C 3 , C 4 , etc). Large molecular 
structures of this kind (>20 carbon atoms) are expected to be flexible enough to form 
an aromatic ring structure. One important difference between this soot formation 
scheme and that involving acetylene and PAHs discussed above is the presence of a 
large number of unfilled valences at the periphery of the planar aromatic structure. 
Therefore, even more than in pyrolysis of hydrocarbons where H can (temporarily) 
saturate these dangling bonds, minimizing this "surface free energy" is expected to be 
very important in soot formation in R Cr B and WC 8-10 stars. This may act to 
enhance the value of pentagons in the aromatic structure (eg., minimize the number of 
surface bonds). The introduction of pentagons induces warping and curling of the 
aromatic planes, which may lead to closure of the aromatic structure on itself, again 
minimizing the number of dangling bonds. Indeed, laser vaporization studies of 
graphitic materials under a H poor atmosphere clearly illustrate the importance of 
such considerations (96,97). 


A second difference between the outflow from these objects and that from C-rich 
giants is the presence of ionizing radiation. WC 8-10 stars are very hot stars and 
most of their energy is emitted in the far ultraviolet. Indeed, the wind from these 
stars is highly ionized. The observed column density of neutral C is typically only 
10' b of the elemental C. R Cr B stars are much cooler (cf., table 3), but they do 
possess an active chromosphere, emitting FUV radiation which may lead to ionization 
(as well as photodissociation). In thjs case simple, neutral, carbon bearing molecules 
do, however, dominate the composition of the envelope. Because ion-molecule 
reactions are much faster than neutral-radical reactions, ionization may play an 
important role in soot formation, in particular in WC 8-10 stars. The bottleneck in 
sopt formation is now the formation of the first molecules. In view of the high degree 
of ionization, the absence of H, and the low pressure, radiative association of C with 
C + and C 2 + are expected to be important reactions in the outflows from WC 8-1 0 
stars- Dissociative recombination of C3 + will lead to Co which can then be used as a 
buildmgblock towards larger acetylenic and aromatic species through ion-molecule 
reactions. Because of its low formation temperature, the carbon dust is expected to 
be highly disordered. Indeed, the IR spectrum of Ve 2-45 shows a weak emission 
Mature at about 7.7 pm, which is attributed to the C-C stretch of aromatic carbon 
(112). This feature is only expected in disordered carbon with an aromatic domain 
size less than about 50A (1). 

F) Carbon stardust formation in supernova ejecta 
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It is appropriate to consider briefly the formation of carbon dust in supernova ejecta 
both because of the (possible) contribution of supernovae to the carbon budget of the 
qalaxy (cf., table 2) and the presence of a carbonaceous component in meteorites 
containing isotopic anomalies characteristic of supernovae (3,1 1 3). Although IR 
excess emission has been detected in several supernova remnants, this may actually 
result from preexisting dust which "lights" up due to the SN (eg., infrared echo) and 
there is presently no conclusive evidence for the formation of (carbon) stardust in bM 
ejecta (47). Because of its timely nature, this discussion will be mainly geared to the 
recent supernova in the Large Magelanic Clouds (SN 1987a). 


Supernova 1 987a in the LMC is thought to result from the explosion (~1 0^ erg) of a 
sliahtly metal-poor star with an approximately 6 M 0 He core and a 5-1 0 M® h 
envelope (114,115). The main-sequence progenitor, SK-202-69, probably nad a mass 
of about 20 M n and had undergone considerable mass-loss, presumably in a proceeding 
red supergianf phase, as evidenced for example by the detection of a nitrogen-rich 
circumstellar shell (116). For our purposes we only need to consider the He core in 
which static He burning has produced a large mass fraction of ^C. Explosive 
nucleosynthesis following the passage of the shock front will convert a large fraction 
of this core into heavier elements, but a carbon-rich zone containing about 2 M 0 
remains (115). The competition between different heat sources (diffusion of 
radioactive heat due to decay of freshly synthesized 56 Co; H recombination; possible 
pulsar) and sinks (adiabatic expansion; radiation) hampers an acurate prediction of 
the temperature-pressure history of the He core material. Uncertainties in the 
opacity due to velocity shearing, compositional gradients, as well as mixing of 
different zones due to Rayleigh-Taylor instabilities further compound this problem 
(115). Indeed, estimates of the gas pressure in supernova ejecta, when the kinetic 
temperature reaches the typical dust condensation temperature of about 1000K, vary 
by aDout 10 orders of magnitude (117). As a result of these uncertainties, predictions 
of grain formation in SN 1987a resemble a weather forecast: it may rain grains or it 
may not. 


G) Summary 

This detailed examination of the chemical pathway towards soot has identified 
several important physical concepts and, although there is some "loose" 
correspondence, tne distinct differences between the kinetical and thermodynamical 
approach should be stressed. Clearly, the role of hydrogen in soot condensation cannot 
be overemphasized. Hydrogen saturates the free valences at the periphery of the 
aromatic nucleation centers and thereby reduces its surface free energy. As a result 
surface free energy plays much less a role in dust formation than envisioned by 
classical nucleation theory. That C 0 H 2 rather than C is the gaseous precursor 
molecule to soot is a further distinction between the chemical and classical 
approach. The saturation with H implies that kinetic effects are very important, since 
a peripheral H has to be abstracted before further chemical growth can occur. Such 
kinetically controlled growth will lead to a highly disordered structure of the 
resulting condensates. The dominant role of H in soot formation is exemplified by the 
high abundance of PAHs as well as the high H content of soot. Physical growth (ie., 
clustering) can also play an important role in determining the final size and gross 
structure of the product (cf., fig. 7). At what point clustering of microcrystals 
overtakes chemical growth by monomer addition will depends on the detailed 
temperature and density history of the cooling gas. Finally, this discussion pertains 
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PROJECTILE 



(a) 



Figure 1 0 . A schematic representation of the impact of a projectile on a target. The fate of the shocked 
target material depends strongly on the peak shock pressure. Since the strength of the shock wave driven 
into the target decays with distance, the impact point is surrounded by successive zones in which the target 
has been vaporized, molten, or converted into a high pressure phase (eg., diamonds). The arrows indicate 
the velocity vectors in the excavation flow which will carve out the final crater (see text for details). 
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to soot formation in a hydrogen-rich atmosphere. When H is absent, as in WC 8-1 0 
stars, the chemical pathway towards dust will be quite different, leading perhaps to 
a different soot structure. Ions or Buckminsterfullerene rather than PAHs may be the 
reactive intermediaries. Nevertheless, it is expected that kinetic rather then 
thermodynamic effects will control soot formation and the structure of the resulting 
grains. 


VII Grains and Shocks in the Interstellar Medium 


Strong shock waves generated by supernova explosions are the dominant destruction 
mechanism for interstellar dust grains (118-121). Behind the shock front, the charged 
grains gyrate at high velocities around the magnetic field lines. Impinging gas atoms 
can sputter atoms from the grain surface (non-thermal sputtering), while collisions 
with other grains can lead to partial and complete vaporization. Evidence for 
destruction of interstellar grains by shock waves is provided by the observation that 
high velocity (v»100 km/sec) interstellar gas (presumably recently shocked) has 
much higher elemental abundances in the gas phase than the general interstellar 
medium (122,123). Although grain destruction by sputtering and vaporization have 
received the greatest attention in the literature, shattering and phase 
transformations due to grain-grain collisions are equally important. The study of 
phase transformations due to grain-grain collisions in interstellar shocks has 
recently gotten additional impetus by the discovery of small («50A) diamond grains in 
carbonaceous meteorites (124). This carbon phase (socalled C-8), which comprises 
about 2% of the elemental carbon, contains isotopic anomalies which presumably 
originated outside the solar system (3). These meteoritic diamond grains may result 
from shock processing of interstellar graphite or amorphous carbon grains. High 
velocity (>8 km s' 1 ) grain-grain collisions behind strong supernova shock fronts are 
expectea to produce the high pressures required to transform these carbonaceous 
materials into diamond. Such collisions are predicted to produce polycrystalline 
diamonds with microcrystal sizes ranging from about 5 to 100A, possibly intermixed 
with some highly disordered diamond (125). 


Behind the shockfront small grains (<500A) are rapidly dragged to a halt through 
collisions with the gas. Essentially, because of their larger volume to surface ratio, 
large grains will not be stopped as quickly (cf., 121). Extinction observations show 
that interstellar grains have a steep size distribution ranging from about 3000A to 
molecular sizes (1 3-1 7). Most of the interstellar grain volume is in the large grains, 
but the surface area is in the small ones. The most important grain-grain collisions 
behind a shock are thus due to a large grain impacting a small grain at high velocity. 

A) Grain-Grain collisions 

Initially, when a projectile impacts a target at high velocity, strong planar shock 
waves are driven into both target and projectile at nearly constant pressure and 
velocity as the kinetic energy of the collision is thermalized. This first stage ends 
when tne shock wave in the projectile has reached the back surface and is reflected 
as a rarefaction wave which eventually reaches the target. These rarefaction waves 
as well as the expansion of the shock wave from the impact site will transform the 
planar shock in tne target into a more or less hemispherical form (eg., the detached 
shock; see figure 10) and its strength will decrease with distance from the impact 
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site. In this second stage the shock pressure will drop approximately inversely with 
the engulfed volume (p-v 1 - 1 ) where the small deviations from the solution for a 
Sedov-Taylor blast wave are due to the backsplash of ejected material which carries 
away energy and imparts momentum (126,127). The rarefaction waves from the free 
surfaces will set up a flow pattern for the shocked material upward and away from 
the impact site (eg., the excavation flow; see figure 10). During the excavation of the 
crater, some of this material reaches the free surface and can be ejected as small, 
shattered fragments. The remainder will line the crater. Shattering may also occur 
when the shock breaks through the backside of the target and is reflected backwards 
as a tensile wave (ie. , spalling). 


The rarefaction wave that follows the shock wave will release the compressed 
material and during this expansion the shocked material may undergo a phase 
transformation. For very strong shocks, the pressure is high enough to cause 
vaporization. For weaker shocks, the material can unload into the liquid phase, or into 
a high density solid phase such as diamond. A shock slows and weakens as it expands 
through a target, so that an initially very strong shock successively vaporizes, melts 
transforms and shatters the target (cf., figure 10). Clearly, the process with the 
smallest threshold pressure will affect the largest volume. Assuming the 
experimentally measured threshold pressure for the graphite-diamond transition of 
600 kbar (i.e., threshold energy e=0.7eV/C atom; cf., § B) and for vaporization of 5.7 
Mbar (i.e., a binding energy e=7.5eV/C atom), the ratio of the volume of diamonds 
produced relative to that vaporized is about 7.5 for semi-infinite targets. Of course, 
for high velocity impacts of equal-sized grains both grains will be vaporized, but 
those collisions are fairly rare. Nevertheless, for interstellar grain-grain collisions 
a careful analysis has to be carried out taking into account the expansion of the shock 
wave in the target as well as the finite-size of the grains. The threshold collision 
velocity for diamond formation is the sum of the particle velocity behind the shock in 
the target and projectile (=2(2e/m c ) 1/2 ; m c =mass of C atom), corresponding to 6.8 km 
s' . Likewise, the threshold collision velocity for vaporization is calculated to be 
about 23 km s' 1 assuming a threshold energy equal to the binding energy. Actually, 
because of kinetic effects, the shock will have to supply several times the binding 
energy to cause appreciable vaporization (127,128) and this estimate as well as that 
of the relative volumes is only a lower limit. 


In order to estimate the effects of grain-grain collisions in interstellar shock waves 
on interstellar grains several different steps have to be taken. First, for a given size 
ratio of projectile and target grain the diamondized volume (or vaporized/shattered) 
is calculated as a function of relative grain-grain collision velocity, taking the finite 
size of target (and projectile) into account. This volume is averaged over impact 
parameter, since glancing collisions couple less of the projectile f s energy and 
momentum to the shock wave in the target (127). This is tnen used as input for a 
numerical code which calculates the structure of an interstellar supernova shock 
wave (120,129). Behind a supernova shock front different sized grains have a 
different velocity history and thus for each grain size at each spatial point in the 
postshock gas the collision rate with all other grains has to be weighted with the 
grain volume affected by such a collision. Integration over position behind the shock 
yields then the diamondized volume for each grain size. The total grain volume that is 
diamondized by a supernova shock at one particular shock velocity is then calculated 
by integrating over grain size. This calculation is repeated for several different 
supernova shock velocities. Finally, these volumes are convolved with the frequency a 
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grain is shocked by a supernova shock with a particular velocity to calculate the 
lifetime of a grain versus diamondization (or vaporization or shattering). This atter 
step requires a model for the interstellar medium. Evaluation Pf these integrals 
x/ipfris an ftxoftcted fraction of elemental carbon in small diamond grains Q about 5 A 
(127) This fraction is not very dependent on the absolute supernova shock frequency 
since supernova shocks both form and destroy (through sputtering and vaporization) 
diamonagrains. Of course, the large number of other steps required in the evaluation 
do make tnis number quite uncertain. Nevertheless, since observations of elemental 
depletions in high velocity gas testify to the destruction of interstel ar grains by 
supernova shock waves, a small fraction of the elemental C (~5 /o) is expected to be in 
the form of interstellar diamonds produced by grain-grain collisions in the 
interstellar medium. 

B) Diamond metamorphism 

Numerous laboratory experiments have shown that graphite and amorphous carbon 
transform into polycrystalline diamond (crystallite size=200A) by a diffusionless 
process when a shock wave traversing these solids drives the pressure above about 
400 kbar. The transformation is complete when the pressure exceeds about 600 kbar 
(130-135). Experiments have been performed on fairly well ordered forms of 
Dvrolytic grapnite (i.e., polycrystalline carbon), on slightly disordered, natural 
graphite from Ceylon, as well as on highly disordered and porous forms of graphite. 

For the more disordered forms the threshold pressure for complete diamond 
transformation is actually somewhat lower (=400kbar). Because of the high 
activation energy (=7.5eV/C atom) for the diamond-graphite rearrangement, 
araphitization of newly formed diamond will not occur whan the cooling behind tha 
shock front is very rapid. Indeed, experiments have shown a high conversion 
efficiency. Even at shock pressures of only 250kbar, weH below the complete 
transformation threshold, about 25% of diamonds have been recovered (135). A single, 
perfect crystal of graphite, when shocked along the c-axis, actually transforms into a 
hexagonal form of diamond (i.a., lonsdaleita) rather than the cubic form. 


diagram : Figure 1 1 shows the high pressure phase diagram for 
ie different regimes can be discerned, where graphite, diamond < 


carbonic regimes can be djscerned, where graphite, diamond and 

the liquid phase respectively are thermodynamically stable. The liquid-diamond 
transition nas been shown with a positive slope, reflecting recent theoretical and 
experimental studies (137,138). Although diamond is thermodynamically favored at 
hiqh pressures and temperatures, grapnite is metastable in part of this regime, 
presumably because the transformation entails a complete rearrangement of the 
carbon lattice. Likewise, diamond is metastable in part of the graphite domain 
Measurements of the metastable graphite melting curve within the diamond stability 
reqime (eg., between about 2500 and 4000 K; 136) have been boldly linearly 
extrapolated to much lower temperatures in figure 1 1 . Theoretical calculations show 
that a metallic carbon phase with a simple cubic structure is preferred at very high 
pressures (>10Mbar; 139). Although preferred at high temperatures and low pressures, 
the exact location of the vapor phase is not well known. Estimates of the 4 

solid-liquid-vapor triple point range from about 0.1 kbar and 4500K to about 2x10 
kbar and 4000K (136). Essentially this uncertainty reflects a debate on the existence 
of a fourth polymorph of carbon which may have been detected in some carbonaceous 
meteorites and in chaotite from the Ries impact crater (140,141 ). In this solid, 
known as carbyne, the carbon is thought to be arranged in linear chains either with 
alternating triple and single bonds (acetylene polymer; cf., figure 2) or with 
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CARBON PHASE DIAGRAM 



g fl ure V - The carbon phase diagram at high pressures and temperatures (adapted from reference 136). 
The regimes in which graphite, diamond and liquid carbon are thermodynamically favored are indicated. 
Graphite is metastable below the graphite melting line within the diamond stability regime. Diamond is 
metastable below about 2000 K in the graphite stability field. The line labeled shock Hugoniot connects the 
loci of shock states accessible for graphite. The stronger the shock the further to the right the immediate 
post shock conditions will be. The diamond symbol indicates the experimentally determined threshold 
conditions required for the graphite -> diamond transformation. 
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cumulated double bonds (cumulene; 142). However, although carbon in the form of 
acetylene polymers could certainly be metastable in a large part of the phase 
diaaram the existence of this fourth stable carbon polymorph has never been shown 
directly Moreover, the original discovery is very much in doubt, since later research 
on several of these samples either discovered impurity inclusions in these materials 
or could not repeat the earlier results (143). 

The oraDhite shock Hugoniot (144) has been superimposed on the carbon phase diagram 
fcf fiaure 1 1 ) The Hugoniot is the loci of shock states accessible for graphite. The 
higher the shock velocity, the higher the resulting ppstshpck temperature and 
pressure will be and the farther to the right on the Hugoniot curve the material 
immediately behind the shockfront will find itself. After the passage of the 
shockfront the compressed material will expand again. This expansion is adiabatic 
and no entropy change will occur. For low postshock pressures, the expansion isotrope 
will almost coincide with the shock Hugoniot. For the high shock pressures 
appropriate for this discussion, the material can still contain appreciable internal 
enerqy after expansion and the expansion isotrope lies to the right of the Hugoniot. 

The approximate location of the shock conditions for which the graphite-diamond 

transition is first observed is indicated by a diamond in i figure 1 

production of (shocked) diamonds occurs typically up to about 1 Mbar and 3000K 

(136,145). 

The onset of the shock generated, graphite-diamond transformation occurs close to 
the extrapolated graphite melting line in the carbon phase diagram and is probably 
initiated by pseuao-melting of graphite and rapid recrystallization in the diamond 
phase (136,146). Extrapolation of these laboratory studies to interstellar grain 
collisions is somewhat nampered by the difference in scale size. Laboratory 
determined crystallization timescales are faster than 10'® s (135), but this may 
still be longer than the adiabatic shock cooling timescale in an interstellar grain 
(Ar/C c *=10‘ s w here Ar (=1 00A) and C s («5 km s'* ) are the size and sound velocity 
in the compressed material) and trie resulting crystal may be highly disordered. Note 
that although the shocked region — ** *h,n tho cmaif ,ntfir«teii?ir nram 

size scales, radiative losses due 1 

important since the radiative lifetimes (for — . . - - ■ - - 

respectively; 147) are much longer than the adiabatic expansion timescale. 

Mir rnr.rvstal sizes : Near the melting line the structure of monatomic liquids is that 
obtained by dense random packing of hard spheres (148,149). Energy and strain 
considerations favor then a local liquid structure consisting of 20 tetrahedra 
arranaed in clusters of icosahedral symmetry (containing 13 carbon atoms; 1 50,151 ). 
These abundant icosahedral clusters can act as nucleation centers for crystallization 
into an ordered crystal structure. The rate at which atoms in the liquid-solid 
interface move to lattice sites is limited by the sound velocity (eg., vibrational 
timescale -10'™ s). The advancing liquid-solid interface will be driven by the 
liberated latent heat of melting at velocities up to the sound velocity (explosive 
crystallization; 152). That is, 

dr 
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where K is the heat conductivity (=3 W cm' 1 K' 1 for diamond at 1 000K; 1 53) and 
1*1 u the temperature gradient (=10 9 K cm' 1 for a 100A compressed region). The 
latent heat of melting of diamond, AH m , is not known. Assuming it is the same as for 
graphite («25 kcal/mol; 136) yields Vj*i km s* 1 . Crystallization will stop when 
neighboring diamond crystals touch or when the heat loss to the surroundings 
becomes too large. The remainder of the carbon liquid will then cool so rapidly that it 
will solidify in a highly disordered diamond structure (c.f., 1). Probably, this will 
become important when the compressed material has expanded again due to the 
rarefaction wave. The diameter of the (ordered) diamond crystals grown this wav will 
a ^ er ^ e J? e1 ky ttj® crystal-liquid interface velocity, Vj, and the time available, (ie. 
Aa=2v/c s Ar). This limits the size of a diamond microcrystal to the size of the 
expansion region behind the shock (Ar), which is typically about 10% of the depth of 
the shock front in the target (1 54). For the formation of interstellar diamonds 
collisions of small (100A) grains with large carbon (a=1000A) grains are most 
important resulting in a Ar of 1 00A. For the parameters quoted above we find that Aa 
is typically 40A. Note that this is to first order independent of the size of the 
expansion region, essentially because the crystal-liquid interface velocity adjusts 
itself to the size of the expansion region. It is therefore not very surprising that the 
sizescale of the microcrystallites in macroscopic experiments (=100A' 135 145 155) 
is very similar to those expected in the collision of submicron particles. 

Once the liquid has solidified, further rearrangement of the lattice will be inhibited 
by a considerable activation barrier. Therefore, the bonds between carbon atoms in 
the interfaces of neighboring microcrystals will be higly strained and distorted. As a 
result these interfacial carbon atoms will show some degree of sp 2 bondina (cf fia 
2), J| IV ® pPP arent "amorphous carbon” character to some types of spectra 

(cf., 124,156). Both laboratory produced diamonds as well as meteoritic diamond 
reveal the presence of such highly disordered interfacial carbon (156 . Upon coolina 
and expansion, the shock synthesized diamond will enter the graphite stability 
regime. Graph itization of the newly formed diamonds will, however, be inhibited bv a 
large energy barrier since all the bonds have to be broken and the lattice rearranged 
Diamond is therefore metastable below 2000K in the graphite stability field. Since a ’ 
large fraction of the carbon is recovered as diamonds in macroscopic shock 
experiments of graphite and since cooling will be even faster, grapnitization is 
expected to be unimportant on the submicron sizescale of interstellar qrains We 
note, however that ion-implantation during further shock processing will lead to an 
approximately 5A thick, amorphous, surface layer in whicn the binding has an 
aromatic character (127). y 

Finally, in those cases where spalling is important (see below), a small fraction of 
the compressed material may be released during the spalling process. In that case 
!u u IZ u S resulting fragments will be determined by the surface free energy of 
smaller 1 pressure carbon llc 1 uid and the resulting solidified droplets may be soniewhat 


Shattering : For the prediction of the structure of interstellar diamonds the effects of 

,*!h?5fo ng must fit 0 be considered. The velocity associated with the excavation flow 
which carves out the crater in the target is much less than the velocity of the 
rarefaction wave and thus the molten material is expected to solidify long before 
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shattering at the free surface can occur. It is important to recognize that the 
fracture yield strength is highly size dependent (127,157 158). For macroscopic 
bodies shattering is due to the initiation, growth and coalescence of mi .{j r oscopic 
cracks, but these play no role at micron size scales where matei rial9 » wiH show a 
ductile (fluid) rather tnan brittle behaviour. The dynamic yield strength of 
polycrystalline diamond is measured to be about 400 kbar (159,160). This is very 
similar to the minimum dynamic pressures required to form diamond by 
shock-metamorphism. Since the pressures associated with the excavation flow are 
much less than those of the shock and since most of the shpck diamond synthesis 
occurs at the lowest pressures (V~P'°- 9 ; cf., § VIIB), shattering during cratei* 
formation is expected to be of little importance. Shattering can also occur when the 
sh^ reaches the backside of the target and is reflected as a tensNe wave (eg^. 
spalling). However, in view of the increased tensile strength of mater ate at small 
sizesca es, this plays much less a role than for macroscopic bodies. Indeed while 
"catastroDhic destruction" of cm-sized glass spheres occurs for pressures of the 
otoer of one kbar (161-163), for submicron-sized glass .spheres the 'critical [pressure 
is of the order of 150 kbar (127). For shock synthesized diamonds this critical 
pressure is expected to be around 400 kbar. Assuming that the typical diamond 
Forming collision occurs between a 100A and a IpOOA grain, the minimum eli sion 
velocity for catastrophic destruction is then calculated to be about 100 km 9 .(J 27 )- 
Since diamond formation in the interstellar medium is dominated by low velocity 
shocks (V=50 km s 125), shattering of newly formed diamonds is expected to be 

relatively unimportant. 

ImDurities- Interstellar graphite or amorphous carbon grains may contain considerable 
amounts^ impurities, either interstitial or substitutional In Particular J^ust 
formed in the outflows from carbon-rich red giants may be heav' hydmgenated 
(H/C=0 1 ■ cf , § VIB). Some impurities, such as N, can be substituted for C m the 
diamond lattice and their presence will not have a major effect, hydrogen (or 
deuterium) atoms are small enough that they can be stored in interstitial diamond 
sites. Natural diamonds contain typically about 0.3% of substitutional nitrogen, which 
in fact gives them their striking colours. About 0.2% of interstitial ® nt . nf 

similar amount of chemically bonded hydrogen is also present ( 164 )- s o ^ a l of 

imDurities will not affect the shock Hugoniot or the graphite-diamond transformation 
very much. Indeed, diamonds have been recovered from shock experiments on glassy 
carbon (hiahly disordered C resulting from heat treatment of furfural resins) 

Sntainina up to 1 0% of oxygen, 3% of hydrogen and 0.2% of nitrogen (1 65). n 
principle, larger amounts of impurities might lead to different resuits^evertheless, 
theoretical and experimental studies of liquid benzene shocked to about 700 kbar 
indicate that some of the carbon is transformed into diamond, although no recovery 
studies have been performed (166). However, we do note that small diamonds have 
been recovered in static pressure/flash heating experiments on severaloraanic 
compounds, including anthracene, naphthalene, chrysene and pyrene (1 67). The 

threshold conditions for diamond formation in these e . x P® r ^ even 

somewhat hiqher than for pure graphite experiments (168). Thus, it is likely that even 
a fraction of impure" amorphous carbon grains such as hydrogenated amorphous 
carbon formed around carbon-rich red giants will transform into diamonds under 
shock conditions. 


In the context of the small diamonds discovered in carbonaceous meteorites, the fate 
of interstititial impurities such as Xe HL is of considerable interest. Pulsed laser 
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annealing Judies of ion-implanted silicon layers show for example that due to 

Iquid ph ? se ’ ,m purities are concentrated near the surface of the 
crystals (169). The measured crystal-liquid interface front velocity in these 

? X fn co? b 2 ut 1 m s ’ is . less than the effective "diffusion" velocity 

(~10 m s .169). However, crystallization will occur much faster than diffusion in 

d^stortf the dfamond structure* ' mpUn " eS are expeC,ed '° be incor P or ^ into (and 
C) Diamonds in the sky 

Since only a small fraction of the elemental C is expected to be in the form of 
d'amond grains due to shock processing, it may be difficult to detect them in the 
^erstellar medium. Small diamond grains will show a pronounced far ultraviolet 
absorption edge and may contribute to the far ultraviolet rise in the observed 
interstellar extinetjon cun/e (125). However, a size distribution of large 

Sil a i 'nH ,v a M 7 °m d i 9 T! Wi " pr i duce a rather smooth extinction curve over 
b ? and U Y ( 17 ^- n J hat case> the Presence of diamond grains will be difficult 
to prove (or negate), since the extinction will be dominated by other more abundant 
dust components (cf., table 1). The density of phonon states of diamond shows several 
4 -5 and 5pm, but these modes are IR inactive for a perfert monovalent 
crystal ancJ thus for large, crystalline diamonds. However, some IR activity is induced 
^he disorder introduced by the presence of impurities, defects and pSularly 
microcrystalhte surfaces (1). Thus, higly disordered, interstellar, diamond qrains y 
resulting from rapid quenching of the liquid (cf., § B), will show weak infrared 
signatures at the peaks of the phonon density ot states. 


The interlocking microcrystals in the polycrystalline diamond qrains will be 
f p epara i e h d by highly disordered thin interface regions. The caFbon Itoms in these 
regions have a large sp 2 character (eg., aromatic structure - cf fia 21 an H thpir ir 
spectrum will resemble that of soot. Tnese polycrystalline conglomerates miaht 
^tnbute to the widespread, near-infrared (12pm) cirrus observed by IRAS ?1 25) 

Ili?i t hii h rfnH?i?/° ntraSt PAH mo,ecule s. they absorb over a large part of the 

visible and UV region of the spectrum without showing distinct structure in the 

extinction curve Moreover, their IR emission spectrum will be very similar to that of 
soot and is therefore likely to resemble that observed in reflection nebulae (8) Third 
because most modes are IR inactive, these grains will radiate more easily in the mid 

« lt ar L pt tl e n ™ atena,s Wlth comparable sizes. A high temperature is also indicated 
by the high Debye temperature of diamonds. Moreover, due to the difference n 
and 'binding, it is perhaps possible to localize the excitation! enemy in a 
small (amorphous carbon) part of the polycrystalline grain, which will temporarily be 
t0 ra nK te t ” rough the near ' and mid-infrared modes (171? Obviousy 
localization effects'^ 1 ^ ^ important t0 stud y the reality of such energy V ’ 


VIII Interstellar Carbon Dust and Carbonaceous Meteorites 

One of the most interesting developments within the field of interstellar du<;t in 
recent years is the realization that some interstellar grains may have been 

X r Eitv d m iTi 7 §i nt £ S 2 nd ' nt ? rp l? netar V du s* 9 partides without totally losing 
their identity (3,1 13,172). Evidence for this rests on the measurement of isotopic 
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anomalies in meteorites and in particular in carbonaceous chrondites. Although the 
meteoritic composition is in a global sense remarkably homogeneous, • ■ _i_ 

non-mass-dependent isotopic anomalies do exist for many elements. These include 
the noble oases the light elements (H,C,N, and O), and tne heavy elements (eg-, Ca, Ti, 
Cr Ni Nd, Sm and others). Although some unusual processes in the solar nebula might 
have Drod’uced non-mass-dependent isotopic fractionation in some elements, it is 
unlikely that they could account for all of them. Moreover, the measured ls0t 2£!S.| o 
anomalies are very characteristic for a presolar origin of the material. For example, 
the Xe and Kr isotopic anomalies associated with one particular meteoritic carbon 
phase carry the signature of s processes in red giants, suggesting the °* 

largely unmodified carbonaceous grains in meteorites whicn were produced in 
rarbon-rich r©d oiant outflows (17,25). Thus, motooritic m8t6nals carry a 
nudeogenic record of the birth site of the dust grain, modified nPta?v° h 

occured in the interstellar medium as well as in the solar nebula or on P' an ®^ r y 
bodies. Consequently, laboratory analysis of meteoritic and ^P ^f. a ^ y nt material 
mav vield detailed microscopic information on interstellar dust, which is not 
obtainable by astrophysical observations. It should, however, be em P f ^ z ?S t ^ a ! t 
meteoritic materials are very heterogeneous and contain dust grains from different 
origins; including star dust from many different birth sites, interstellar medium 
dust and solar nebula condensates. Generally, the actual presolar carrier of the 
measured isotopic anomalies is not known. Such information is, however, of prime 
importance for our reading of this record and for our assement of the implications for 
interstellar dust and its evolution. 


Manv different carbon phases with isotopically anomalous composition have been 
discovered in carbonaceous meteorites. Table 6 provides a summary of some relevant 
data on these phases (cf . , 113). Each of these will be discussed in turn. 


Thft C-a Dhase: This chase consists of turbostratic carbon, that is a carbon solid ^ 
containina d i so rienteo graphitic layers quite similar to soot. It contains almost pure 
^Ne. Most of the^Ne in the galaxy is produced by WC stars. These stars also produce 
minor amounts of carbon staraust (cf., table 2) and thus might be the origin of f this 
SrtSJn ’phase in carbonaceous meteorites (1 73). However it is difficult to produce 

nudwjsyntheslstn S^out 0 02(44). An 

^Na), as well as other elements (^g, ^Mg,^A|L 27 AI) must have been rap^gly 
incorporated into the carbon staraust before the 22 Na decayed (t«2.5 yrs) to Ne. We 
note that electron donors such as alkali, alkaline earth and transitio n metals , read ily 
form intercalation compounds with graphitic and turbostratic carbon (175). Unlike the 
noble gases these metals are, however, difficult to measure in meteorites. 


This meteoritic carbon phase component is also known to be arched in J 3 C and 15 N 
(cf table 6) However, the 13 C enrichment is less than one might expect from 
theoretical models for novae nucleosynthesis which predict 12 C/ 13 C ratios 'ass than ^ 
unitv (30 53). Also 13 N overabundances of >10 3 are predicted by some models. These 
predictions are quite sensitive to the details in the models than 

shell burning) and the problem may reside in the nucleosynthesis models rather than 
the identification of this meteoritic carbon phase with novae stardust. Alternatively, 
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Table 6: Interrelationship of interstellar and solar system carbonaceous materials 


Carbon Phase 

Isotopic 

Anomaly 

concentration 

[ppm] 

Birth Site 

C-a 

(Ne-E) 

pure 22 Ne 
' 2 C/ 13 C=70 

5 

-Novae 

c-p 

s process Xe 

l2 C/ 13 C=40 

86 Kr 

<5 

-C-rich 
Red Giants 

c-s 

(diamond) 

Xe-HL 

14 N/ 15 N > 400 
12 C/ 13 C=90 

400 

-graphite -> diamond 
in ISM shock waves 

C-e 

22 Ne 

12 C/ 13 C<10 

2 

-Novae 

soluble and 
dust 

D/H up to 1 0' 3 

20,000 

-pyrolysis of ISM 

insoluble 
organic matter 

mundane C and N 


-PAHs (?) 


the C-a phase may be contaminated with solar system carbon with a normal 12 C/ 13 C 


I he g-g phase : This carbon phase in carbonaceous meteorites consists of SiC qrains 
which contain enriched s-process isotopes of Xe and Kr (cf., table 6). It is also 

Ilf generally accepted that this carbon phase represents stardust 
222*^0^ th ® outfl ° w from carbon-rich red aiants (3,1 13,172), where we note 

Irf* Shi! VSffl f/^® ns f |C f° r SiC, is often seen in such environments 

(cf., table 1 and 2). The 14 N/ 15 N ratio in this carbon phase is solar, which is 

Sl nce the seme dredge-up episodes that enhance also will 
produce 4 N (1 76) Because of its refractory nature SiC is more readily detected in 
stepwise heating of carbonaceous meteorites than carbon dust. The latter is 
however, more abundant in circumstellar and interstellar dust and a "pure" carbon 
similar isotopic anomalies is expected to be present in carbonaceous 
meteorites at about the 50ppm level. 

o a i £? arbonac ®° u s meteorites consists of small diamond 

J?i ns JP 24,15 $’ ? oth the heav y ar| d light isotopes of Xe are enriched 

t0 th f sun > jmplying a presolar origin (3). It has been suggested that these 
I s ™? C ?" omalies reflect the effects of mixing of two different nucfeosynthetic 
zones within one supernova. This may have occured either before condensation or 
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after condensation in one zone by ion implantation of the products of the second zone. 
Avariant of the latter is ion-implantation in grains 

ohase bv the overtaking high velocity, supernova ejecta (177). Recently, it has oeen 
suaaested that these smalf diamond* grains are the result of shock metamorphism of 
carbon stardust in the interstellar medium (1 25 . In that model, the measured 
i^lto^c anomalies reflect the sum of the nucleosynthetic processes j|king i place in 
the various birth sites of the original carbon stardust. Both the small N and the 
absence of a enhancement are difficult to understand in terms °fo single . 
rircumstellar source (eg., supernovae or their red supergiant progenitors) of this 
meteorttic carbon phase but might conceivably result from averaging over many 
different carbon stardust sources. In the latter case, one ^©y vever f x P e 9* , . 
S5K? variations on an individual interstellar grain basis M.OOOA scale) ft should 
be ernohasized that the measured abundances of the isotopically anonnalous Xe are 
measured on uo samples and reflect thus an average over about 1 0' ' d, amend 
microcrystals or about 10 8 typical interstellar dust grains (-1000A). 


Tho Qnrfare lavers of this meteoritic carbon phase seem to contain noble gases 

K£S&"J3SS» !B8w^<W?mSSM » «»« ;» 

the^mass^ the 'element under consideration Note that the drift velocitymredgiant 
outftows is generally very small and ion implantation should be ohimportaot Finally, 
anoreciable drift velocities will also occur in interstellar shock waves and the 
resulting ion implantation will tend to dilute the nucleogemc record from the grai 

birth site. 


Th fi r ^ ohase- The C-€ phase consists of SiC grains carrying an almost pure 2 *Ne 
contaminated 

nmanir matter Two comDonents are lumped together under this nomer (3,1 1 3,172^^ 
Firsf n tffe^cid^so ubte fraction (=25%) contains such organic compounds as amino and 
*cids S^nd the acid-insoluble fraction (-75%) consishng of aromatic 
mnitie*; connected bv short aliphatic bridges and sidechains (eg., kerogen). Both 
Dhases S h^y daterated Such chemical fractionation wilT only occur at very 

Although the cooling solar nebula wi I have experienced 
«?nrh low temDeratures the timescale to reach this equilibrium is far too long to 
obtai na s^z^aWef ractfo natio n (3). High deuteration effects have been observed in 
many gas phase molecules in interstellar molecular clouds and 
orfginate from ion-molecule reactions (178). Large 

fiynacted for iev arain mantles accreted in molecular clouds (179). I nese aeuieranon 
effects result from non-equilibrium chemistry driven by cosmic ray ionization and 
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equilibrium is readily established at interstellar timescales (-10 7 yrs). For this 
f® 3 ® 0 " the hypothesis that interstellar, as distinct from circumsteflar material has 
been incorporated into carbonaceous meteorites without major modifications has 
Sfhni^'SS/ ♦h CCepted t 3, 172 - 18 ?)- The actual deuterated carrier(s) into the solar 
hSion W th?5^ ?? s ph ?® e °f S iS ld u tate - and the importance of the distinction 

! h ® solu P le and insoluble phases is unknown. Possibly, the kerogen results 
whlTrf fh S i e DAS of sma| l '^rstellar PAH molecules and amorphous carbon grains 
where the PAHs carried the deuteration effects due to photochemical processes in the 
interste'lar medium ( 81). Alternatively, it may result from prolonged pySfy^s of 
(deuterated) interstellar molecules either in the solar nebula (3) or in the 
interstellar medium (2). In the latter model the acid-soluble phase mav reoresent an 
™X d '^ st 2 p , in K th9 . UV Photolysis and transient heatingof inte7steltar P grafn 
£!?™L 0S ' P deec ?’ la b°ratory studies show that carboxylic acids are readily formed 

exmlSn MV raH-fv ices 5?’’ H2 ° ?u d C P) by such Processes (182,183). Prolonged 
exposure to UV radiation and transient heating is expected to result eventuallv in the 

a most comp ete removal of elements other than C (and possibly H) aXo fead to a 
structure similar to that of meteoritic kerogen. y ' a 


IX Summary 


conclude tba J interstellar dust is a collection of many different 
2S i ere many different stardust birth sites each with distinct physical 
nf^nhia^r 1103 fopd't'ons a nd elemental and isotopic composition. Many different types 
of obj^ts contribute to the carbon budget of the galaxy, including WC stars 

fhoc^e va ? an £ ®arbon-rich red giants. Although the relative contribution of each of 
these is not well known, it is likely that both low mass and high mass stars 

originates in b wr fl’Vnl''. As, ? n0l ™ ca J evidence shows that carbonaceous stardust 
mfi 9 t«nri?r stars - C-rich r ®d giants, and possibly supernovae. Studies of 

amorffi l tr i? 9ly s . u 99®st that some novae produce SiC and 

balance Ahn nt^n"/ 9 S Si p , robably farbon-rich 9'ants dominate the carbon stardust 
~ a )? n ®. e - A Pout 20 /o of the elemental C is injected in the form of dust arains Sinre 

du t ^nr2inc b th 0 e n/atlC i ns ^" est that about half of the elemental C is locked up \n* 

^5 fiTiphes th© presence of an efficient formation mechanism for 
carborocoous dust in the interstellar medium. One attractive mechanism is UV 

mantiotf °i accreted *cy grain mantles, which may result in organic refractory arain 
mantles (polymeric carbon) or perhaps even amorphous carbon grain mS 9 

stars^nd'wr^ta^cho 8 ir l the stardust condensation zone around C-rich giants, R Cr B 
i S ^ ar ®. have been reviewed. These objects span a wide ranae in Dhvqirai 
^m.^ 0 ? 1103 ^ 0 ^ 1 * 1008 ’ yet car bon stardust seems always to form oncethe P ^ 
temperature drops to about 1 000K. The typical densities inferred for the carbon 

This win fnrm^ 3 * hi USt * as t0 “ 

nSSSiJkli? w Li fo . rr H a P robl ® m for any nucleation theory. Stardust formation has 

classi^n 6 ^ thermochemistry (ie., condensation sequence) and 

® a .j ?*? L nuC u ! on J b ® or y- Tne problems with this approach have been oointed out 
to dust emphas,zed that an y theor y should start by identifying the chemical pathway 


The chemical processes that convert C 2 H 2 into carbon stardust in the outflow from 
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rarbon-rich giants are probably very similar to those occuring during the gas phase 
ovrolvsis of hydrocarbon molecules in flames. PAHs may play an important role in 
this Among the PAHs expected to be abundantly produced in such outflows are 
naphthalene, acenaphthalene, pyrene and coronene. Their detection in circumstellar 
shells will form an important confirmation of these chemical schemes as well as 
suDDortina the ubiquitous presence of PAHs in the interstellar medium. Photqchemical 
Drocessina as welfas shock wave destruction will further influence the chemical 
composition of the interstellar PAH family. The carbon stardust formed in these 
outflows is expected to be very similar to soot and to consist of small aromatic 
moities (ie PAHs) clustered in larger units and connected by aliphatic bridges and 
side chains (ie., tetrahedral carbon). The hydrogen content of these platelets may be 
verv hiah (=10%). These carbon platelets will be clustered in the form of larger 
spherical particles. The basic building blocks of carbon stardust wil thus have a 
structure quite similar to that of hydrogenated amorphous carbon films (HAC). 

The carbon stardust formation mechanism in WC 8-10 stars is not understood very 
well It is likelv that it starts with the formation of small acetylenic species such as 
Co Co and their ions. Indeed, ionization may play a dominant role. An important 
difference with carbon stardust formation around carbon-rich red giants is the 
absence oHiydrogen and hydrogenated molecules in the outflow, Clearly PAHs and 
HAC will not form a good model for this stardust. Actually, Buckminsterfullerene 
fCcnl and its larger analogs, formed by laser vaporization of graphite targets in a 
hyirogen-poor atmosphere, may play an important role in the condensation processes 

around these objects. 

Hiah velocity (v>7 km s' 1 ) grain-grain collisions drive strong shock in carbon 
stardust which provides the high pressures required to transform some of it into 
diamond. This supports the suggested interstellar origin for the recently discovered 
small diamonds in carbonaceous meteorites by providing a feasible interstellar 
formation mechanism. Shock metamorphism of carbon stardust is expected to produce 
Dolycrvstaliine diamond grains. The microcrystalline sizes are calculated to be about 
40A in aood agreement with the measured ones in meteoritic carbon and shock 
experiments. Since the diamond formation is thought to go through a liquid precursor 
phase the structure of the original carbon stardust grains is of little importance. 

Indeed highly ordered graphite, amorphous carbon as well as glassy carbon containing 
substantial amounts of impurities will transform into diamonds under shock 
conditions. 


A wide variety of isotopic anomalies have been detected in different carbon phases in 
carbonaceous meteorites. These carbon phases are thought to consist of interstellar 
or circumstellar dust grains which have been incorporated into the meteorite without 
major modifications. The isotopic anomalies reflect then the particular 
(nucleosynthetic) processes which played a role in the births sites of these grains. 


^Na (“Ne)L ion impiamauon ^wg., aw ni_, w pucu>wj, emu h 

the ISM (eg., D/H). Because many different sources contribute, one of the 
characteristics of stardust preserved in carbonaceous meteorites is a widely varying 
isotopic composition on a submicron scale. In particular, the ^C/^C ratio may vary 
dramatically on this scale depending on the source of the stardust ranging from about 
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zero for WC 8-10 and supernovae, solar (*0.01) for optically visible C-rich Miras 0 3 
for some C-nch giants, to >1 for novae. In the protosofar nebula, stardust from ’ 
d'fferent birthsites may have represented distinct, isotopically enriched reservoirs 
because of their different physical and chemical properties, may have been 
released at different times or in different parts of the nebula. Another signature of 
interstellar dust, which may help recognize presolar dust in meteorites or IDPs is 
shock metamorphism on a submicron scale. ’ 
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